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Electronic structure calculations for layered zincblende semiconductors are described within a 
restricted basis formalism which naturally and non-perturbatively accomodates both crystalline in- 
version asymmetry and cubic anisotropy. These calculations are applied to calculate the electron 
spin decoherence times Ti and T2 due to precessional decoherence in quantum wells. Distinctly 
different dependences of spin coherence times on mobility, quantization energy, and temperature 
are found from perturbative calculations. Quantitative agreement between these calculations and 
experiments is found for GaAs/AlGaAs, InGaAs/InP, and GaSb/AlSb (OOl)-grown quantum wells. 
The electron spin coherence times for CdZnSe/ZnSe II- VI quantum wells are calculated, and cal- 
culations of InGaAs/GaAs quantum wells appropiate for comparison with spin-LED structures are 
also presented. 

PACS numbers: 78.20.Ls, 72.25.Rb, 73.21. Cd, 71.15.-m 



I. INTRODUCTION 

In recent years there has been considerable interest 
in exploiting not only charge but also spin in solid- 
state electronics, which has led to new ultrafast op- 
tical studies of electron spin dynamics in bulk and 
quantum-well semiconductors and possible applications 
to ultrafast spin-dependent switching, quantum com- 
puting, or other "spintronic" devicesi*^*^. The magni- 
tude and persistence of such effects are governed partly 
by the spin coherence times Ti and T2, describing the 
decay of longitudinal and transverse spin order, re- 
spectively. Ultrafast optical measurements have been 
performed of both Ti and T2, although in different 

j,^-^j^^j,j.34.5.6.7.8.9.10.11. 12.13.14.15.16.17.18.19 

We recently described an approach to calculate these 
times for layered zincblende semiconductor structures 
based on fourteen-band non-perturbative restricted basis 
set calculationsiiSiSfliSi. Success was reported in resolv- 
ing the discrepencies between theoretical and experimen- 
tal results for several material systems, both in trends 
and in absolute magnitude. Here we provide full details 
of the calculational method, including both the general 
restricted basis calculations, which can be applied to a 
variety of optoelectronic properties, and the specific cal- 
culations of spin coherence times. 

The mechanism of electron spin decoherence we con- 
sider occurs via the spin precession of carriers with fi- 
nite crystal momentum k in the effective momentum- 
dependent crystal magnetic field of an inversion- 
asymmetric material. This effective k-dependent crystal 
magnetic field is a direct result of the spin splitting of 
the conduction band at finite momentum caused by lack 
of a spatial inversion symmetry in the zincblende lattice. 
Thus a proper calculation of this spin splitting requires 
a proper electronic structure. 

In this work, we employ a non-perturbative het- 
erostructure theory based on a generalized K-p enve- 



lope function theorji2SiSi2iiS^ solved in a fourteen-band 
restricted basis set. The electronic structure of quantum 
wells is obtained by expressing the electronic states as 
spatially dependent linear combinations of the fourteen 
states in the basis. The full Hamiltonian is projected 
onto this restricted basis set, which produces a set of 
fourteen coupled differential equations for the spatially 
dependent coefficient of the basis states (generalized en- 
velope functions). These equations are then solved in 
Fourier space in a similar method to that of Winkler and 
RossleriSSiSi. The expressions for electron spin coherence 
times Ti and T2 are obtained by solving the momentum- 
dependent spin density matrix equation, and the calcu- 
lations of electron spin coherence times are carried out 
within the framework of the fourteen-band theory. 



Our paper is organized in the following way: in Sec.HTI 
a detailed description of the theory of electron spin de- 
coherence in bulk and quantum-well zincblende semi- 
conductors is given. This includes a description of 
the D'yakonov-Perer(DP) formulae for bulk^*'^^''^" and 
the D'yakonov-Kachorovskii (DK) formulae for quantum 
wells^ to facilate the comparison between our results and 
those of DP and DK. We also describe the relationship 
between the longitudinal spin relaxation time Ti and the 
transverse spin relaxation time T2, and the connection 
between the longitudinal spin relaxation time Ti and the 
spin-flip time r^ . The calculation of electronic structure 
and the momentum-dependent effective magnetic field 
based on a fourteen-band non-perturbative K • p theory 
is described in Sec. IIIII Calculated electron spin coher- 
ence times for a variety of material system, both bulk and 
quantum-well semiconductors, are presented in Sec. IIVI 
and comparison with experiments is made where possi- 
ble. 



II. ELECTRON SPIN DYNAMICS 

Spin relaxation of conduction electrons in inversion 
asymmetric bulk semiconductors due to rotation of spins 
during scattering of electrons was originally studied by 
D'yakonov and Pereli22i2S using a momentum-dependent 
spin density matrix approach based on a simple elec- 
tronic band structure approximation (see also Ref . |30|) . 
This work was later extended by D'yakonov and Ka- 
chorovskii^- to study electron spin relaxation in inver- 
sion asymmetric quan tum- well semiconductors (see also 
Refs. [33, I33, and|3J). In this work, we generalize the 
analysis of electron spin relaxation for both bulk and 
quantum-well semiconductors to include multiband ef- 
fects. 



p{t) in the |n,cr, k) representation in which the Hamilto- 
nian H_^ is diagonal: 



ij"|n,a,k}=i?l(k)|n,a,k), 



(4) 



where n and a are the band and spin indices respectively. 
The associated eigenfunction is given by 



(r|n,cr,k) = exp(ik- r)(r|M„o-k}, 



(5) 



where |u„crk} is expressed in terms of the product of the 
Bloch states |m°) at k = and the spin-1/2 states \Ua)- 



Tk) = ^ a„^„/^/(k)|M°,} |C/a'). 



(6) 



A. DENSITY MATRIX FORMALISM 

To study electron spin decoherence in bulk and layered 
zincblende semiconductors, we consider a population of 
Nc electrons in the system. Each electron is described by 
the total Hamiltonian: 



ijT ^H^ + hSO ^j^SC 



(1) 



where H'^ is the periodic crystal Hamiltonian without 
spin-orbit interaction, H^'~' is the Hamiltonian for the 
spin-orbit interaction, and H is the spin-independent 
Hamiltonian describing any scattering interaction. The 
crystal Hamiltonian is H = H° + H^° with H" and H^° 
written as 



In the fourteen-band restricted basis set, u^ G 
{5,X^,y^,Z^,X^r^Z'=} andUa e {T,i} (see Sec. Hn 
for detailed dissussion) . The eigenstates of the electrons 
liy, t) under the action of the total Hamiltonian H^ can be 
expanded in this complete set of time-independent states 
|n,cr,k): 



,t) = ^c^„^(k, i)|?i, CT, k). 



(7) 



naW 



and the matrix elements of the time-dependent spin den- 
sity operator p{t) are then given by 

Pn'a'-k'nakit) = (n', Cr', k' |p(t) |n, CT, k) 

= ^E<"'-'(k',^)c.n.(k,i). (8) 



H° 



2?77,„ 



i/so^ 



4m^c' 



Vir), 



■ [vy(f)xp].^. 



(2a) 



(2b) 



where rric is the free electron mass, c is the velocity of 
light, h is the Planck's constant divided by 27r, p is the 
electron momentum operator, ^(f ) is the periodic crystal 
potential, and & is the Pauli spin operator. We consider 
the case where the spin-independent scattering potential 
is weak (iJ^o ^ ^jsc ^ ^oj g^ ^1^^^^ ^oth ff^O and ff^C 

can be treated as small perturbations. To describe the 
electron spin dynamics in bulk and layered zincblende 
semiconductors, we start from the von Neumann equa- 
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H^,m 



(3) 



where p{t) is the time-dependent density operator for the 
electron population. 

To derive the spin-dependent distribution function for 
the electrons, we shall express the spin density matrix 



The spin density matrix equation for Eq. ||21l in this 
representation is 

d r - 1 

i^-^Pn'a'k'nakit) ^ {n\a',k'\ H'^ , p{t) \n,(T,k). (9) 

Substituting Eq. ^ into the right-hand side of ^ and 
expressing each operator in terms of individual matrix 
elements we have 

d_ 
ot 



= 2^ iEj^,^,(k')dn'n"Sa'(T"S]i'k"Pn"a"k"na]iit) 
n" a"\<i" 

— Pn'cr'k'n"CT"k"(^)-^no-('*)^n"«'^o-"cr^k"kl 

+ 2_^ I ^„'o-'k'n"(T"k"P""o""k"Ti(Tk(fc) 
n" a"]<i" 

+ 2^ \Hn'cr'Wn"a"wPn"(y"'k"na\s.{t) 



n"fT"k" 



Pn'o-'k'n"cr"k"(,ij-W„//o-"k"nCTk 



(10) 



For compactness, we may rewrite Eq. (|10ll in matrix form: 



to is given by 

Pn'k'nkito = -OO) 



{n',k'}7t{n,k} 



= 0. 



(14) 



7 , yM.n'k'n"]i"Pn"k>'nk[t) Pn'k'n"k" (OSn"k"ri 
i"k" ~ ~ 

/ , (^n'k'n"k"Pn"k"Tik(0 ~ Pn'k'n"k" (0^n"k"rik ) J 



n"k 



ra"k" 



(11) 



To study the electron spin decoherence we need to evalu- 
ate the spin-dependent distribution function for the elec- 
trons, which is related to the diagonal elements of the 
spin density matrix Eq. H12|l with respect to indices n 
and k. Inserting Eq. H13I) into the right-hand side of 
Eq. I|12|l for the diagonal elements we obtain 



where [£^° (k)] „ = E'l^,,(k)Sa"a arc the matrix ele- 
ments of E_^{'k), and Pn'k'n'k{t) is a 2 X 2 matrix with 

respect to the indices a' and a. It is convenient to sepa- 
rate the diagonal terms with respect to the indices {n, k} 
and {n',k'} from the off-diagonal terms in the summa- 
tions of Eq. Ullfl and rewrite it as 



ih 



-KTPn'-k'nkit) = (|„'(k') - £„ (k) j p„/k'„k(0 
-[Kjk'n'k'gn'k'nkit) ~ £«'k'nk(i)^nk«k j 
-{M.n'k'nk£jikni<.{t) ~ gn'k'n'k' {t)M.n'k'nUJ 

yKn'k.'n"k"gn"i<."nk{t) ~ gn' k' n" k" {t)Kn" 



n"k' 



SO-SC 
k"nk 



(12) 



where |„(k) ^ ^O(k) + ^^^k, Kl'^k^^ ^ ^fX'„k + 
iin'k'nk' ^^'^ ^^^ prime on the summation denotes the 
sum is over {n", k"} / {n, k}, {n', k'}. Note that K^^^^ 
commutes with Pn'k'nk{t) because ii^^nk i^ ^P^"^ inde- 
pendent and hence it is also diagonal with respect to the 
indices a' and a. In addition, the matrix e;„(k) has identi- 
cal diagonal elements in the |n, cr, k) representation (i.e., 
the spin- up and spin-down eigenstates are degenerate), 
therefore the first term on the right-hand side of Eq. H12() 
can be rewritten as {Sn'a'O^') - £na{^)) Pn'k'nk{t). 

Because the scattering potential H^'^ is weak, the di- 
agonal elements of the density matrix are larger than 
the off-diagonal elements by a factor of the scattering 
strength (see Ref. 1361 for detailed discussion). Further, 
the spin-orbit scattering H^'~' is usually much weaker 
than the momentum scattering H^'~^ in zincblende semi- 
conductors; therefore, we may neglect the second and the 
fourth terms on the right-hand side of Eq. 1)12(1 and ob- 
tain the solution for the off-diagonal elements for the spin 
matrix Pn'k'7ik{t): 



Pn'k'nk{t) 



{71' ,k'}=i{n,k} 



"'' [Jj-n' k'7ikP7iknk[T^ ) — Pn'k'ri'k' (,1 Jiin'k'nk 



X exp 



^(£„.,,(k')-e„,(k))(i-t' 



(13) 



where we have assumed that the initial condition at t 



d 



dt^" 



jrP7ik7ik{t) = — ^„knkiP"knk(i) + —J y ^ 



-1 



n'k' 



dt' 



M:nk7i'k' yU,7i'k'7ik£n^nk[t ) - £„'k'n'k' (.1 j^„' 



SO-SC 
k'nk 



X exp 



— [e,,,„,{k')-ena{k))[t-t' 



,/rrSO-SC„ /,'\ „ /./\ rrSO-SC \ rrSO-SC 

-r|i3.nkri'k'Pn'k'n'k'(,t ) — p7iknk[T^ Jiinkn'k' hS-7i'k'nk 



X exp 



^(e„v'(k')-£«a(k))(i-t' 



(15) 



Taking the limit ?i — » 0, we obtain the kinetic equation 
for the diagonal spin density matrix 

-p„k(i) = ^tek,Pnk(i) 



27r 



-sc 

ik 



E' rrrSO-SC ^ (,\ttSO-SC 

\M.7ikn'k'P7i'k' [T^)M.7i'k'7i 

n'k' 

1 f ttSO-SC ttSO-SC ^ (i-\ I r, {i-\ ffSO-SC rrSO-t 
~2 \^=nkn'k'^n'k'nk£nk(,';j + £nk l,rj^„kn'k'^n'k' 

x,5(£„-,.(k')-£„.(k)), (16) 



-sc 

:'nk 



where p„k(i) = P7iknk{t)- In deriving Eq. (|16|l we have 
made use of an identity for the delta function 



9(—t) lim — exp 



— {E + ie)t 



TrS{t)S{E) + iV 



1 

E 
(17) 



where 9{t) and 6{t) are the usual step and delta funcitons 
respectively, and V[l/E] denotes the principal value of 
1/E. 

Multiplying both sides of Eq. (|16|l by the Pauli spin 
matrices a^ and taking the trace with respect to the spin 
indices a G {Tii}: we obtain the momentum-dependent 
spin-polarization vector 



l="-'" = T^* 



-xE- 



n'k' 



K7ikig7ik{t) : 

rrSO-SC ^ (,\ ttSO-SC 

ilnkri'k'Pn'k' [^ ) IL71' k' 7ik 



1 

~2 
xS 



f ttSO-SC rrSO-SC „ (A 1 r, {i-\ f/SO-SC rrSO-SC 
1 ^nkn'k'^n'k'nkPnkl.rj + Pnk l.tj^nkn'k'^n'k'nk 



(18) 



where Tr stands for trace and the momentum- 
dependent spin-polarization vector is defined as s„k(i) = 
Tr[p„k(0^]- This general kinetic equation is now further 
approximated depending on the relevant situation. 



B. SPIN-INDEPENDENT SCATTERING 

Electron spin relaxation in zincblende-type semicon- 
ductors near room temperature is dominated by the 
D'yakonov-Perel' (DP) precessional mechanism22iS2i^, 
which is a direct result of the spin splitting of the con- 
duction band that occurs at zero external magnetic field 
due to spin-orbit interaction in inversion asymmetric 
crystals. This can be viewed as having an effective 
momentum-dependent magnetic field H„(k) acting on 
the electron spins, which causes Zeeman splitting be- 
tween the spin states Eni{]i.) and £'„|(k). Consequently, 
non-equilibrium electron spins at finite momentum k will 
precess about these effective momentum-dependent mag- 
netic fields at the Larmor frequencies proportional to the 
spin splittings. The presence of scattering causes random 
fluctuation of spin precession frequencies and directions 
leading to decay of the non-equilibrium electron spin pop- 
ulation. 

To study the electron spin dynamics in zincblende- 
type semiconductors, we shall consider only the dominant 
scattering mechanism (spin-independent scattering) in 
Eq. I|16|) pertaining to the DP spin relaxation mechanism. 
For spin-independent scattering (i.e., the off-diagonal cl- 
vanish), equation fH^ reduces to 



ementsofi£SO,„k 



l^-W = T^' 



Knk^Pnk{t) 



- E Tr[E™.(k, k') (^„k(i) - £„'k'(t)) ^] (19) 

n'k' 

where lil„„/ (k, k') is the transition probability of electron 
scattering and is given by 

^(f2'(k,k') =^„„'(k,k')^(e„v'(k')-e„.(k)), (20) 



andWZ{kX)^mir,,.'H^^' 



'k'nk* 



C. KRAMERS DEGENERACY 

Time reversal invariance of the crystal Hamiltonian H 
requires that a spin state at a given k must have the 
same energy (Kramers degeneracy) as that of the op- 
posite spin at — k: £'„f(k) = £'„j^(— k). For spatially 
inversion-symmetric crystals (diamond crystal structure) 
such as Si and Ge, the states of the opposite spin at k 
are degenerate En^ (k) = Eni (k) . On the other hand, for 
spatially inversion asymmetric crystals (zincblende crys- 
tal structure) such as GaAs and InAs, -E„|(k) 7^ _E„|(k). 
Consequently, there is a nonzero spin splitting [Aii'„(k) 



= £'„|(k) — EniCk)] at finite momentum k due to spin- 
orbit interaction. The momentum-dependent effective 
spin-splitting matrix Hamiltonian due to any spin-orbit 
coupling can be written as 



m:^ ^ -n{n,k) ■ ^, 



(21) 



where the momentum-dependent vector ri(n,k) is inter- 
preted as the precession vector whose length equals to the 
spin precession frequency (Larmor frequency) and direc- 
tion defines the spin precession axis. One consequence of 
Kramers degeneracy is that the precession vector J7(n, k) 
is an odd function of k: 



J7(n,k) = -rj(n, -k) 



(22) 



It is obvious from Eqs. H21() and (|22|l that the spin-orbit 
matrix Hamiltonian Hf,^_ is also an odd function of k. 



D. ELECTRON SPIN DENSITY 

In zincblende-type semiconductors, the spin coherence 
time for thermalized electrons is typically much longer 
than the momentum relaxation time; therefore, it usually 
satisfies the condition 



T(n,k)f7(n,k) < 1, 



(23) 



where T(n, k) is the momentum scattering time (see 
Sec. Ill F|l . This allows us to solve Eq. ((T^ by mak- 
ing successive approximations to the spin density matrix 
Pn\s.{t) in the small parameter A(n, k) = T(n,k)ri(ri, k). 
Expanding the spin density matrix Pn\^(t) in powers of 

A(ri, k), the power series of the density matrix may be 
written in the usual form as 



M( 



.(1) 



.(2) 



(m). 



Note that A(n, k) ex H^^ [see Eq. ^^\. It was shown by 
D'yakonov and Perel'iSSiSa (see also Refs.lHandlsl) that 
it is sufficient to expand the spin density matrix to first 
order in A(n, k) (i.e., m = 1). After substituting Eq. l|^ 
in H19|) . we have 



d_ 
di 



Snk(0 



h 

E 



E 



-Tr 



rrSO (0).,N 
^nk 1 £nk y^) 

dk' 



-Tr 



^so „(i) 



S-nk ' 



f^n^it) 



(27r)3 
dk' 

(M 



Tr 



rTr 



Eitl^(k,k')(^i°^(t)-^a(i) 



(25) 



(0), 



where the zeroth order term p^ (t) depends only on en- 
ergy £„CT (k) , while the first order term p^j^ (i) ex H_f^^ . In 



Eq. H25() we have converted the summation over k' into 
an integral ^j^ -^ {2tt)~^ J dk. In the first order ap- 
proximation, the spin density matrix can be written as 



Pnk{t) = (/Onk(i)) 



k ' E-. 



llnlitl 



(26a) 



matrix p'^^ it) is governed by the spin-spHtting matrix 

Hamiltonian H^, while the relaxation rate of /0„i^(i) is 

dominated by the momentum scattering rate W_nn' C*' ^') 
under the condition (|23|l [i.e., T(ri, k) ^ ri~^(n, k)]. 
Therefore, on the time scale of electron spin relaxation, 
which is of order 17^^ (n, k), Eq. H29|) reduces to2^ 



(pi\^W)k-o, 



(26b) 



where (pnk(i))k denotes the averaging of Pnk{t) over all 
directions of the crystal momentum k. Equation Ij26b|) is 
zero because p^^it) is an odd function of k [see Eqs. H21|l 

and (|22() ]. Note that since p„k(^) depends only on energy, 

it is independent of the direction of k (i.e., pl^y.{t) = 

(£ik(^))k = (£nk(i))k)- Integrating Eq. ^ over the 
momentum variable k and summing over n, we obtain 
the kinetic equation for the electron spin density S{t): 



l«w^> 



E 



dk 



(27r)3 



TirSO A^)(j-\ 
U,nl,,PnkW 



(27) 



where the electron spin density is defined as 
rfk 



m 



E 



(2^)^ 



rS„k(0 



Tr 



= Tr 



E 

n 

E 



dk 



(27r)3 \!= 



f^lit) 



d2it) 



dk 



(0)/ 



J^tJm 



(28) 



Equation (|27|l is obtained by noting (i) the averaging of 
the first term on the right-hand side of Eq. (|25l) over 
momentum k vanishes because the spin-splitting matrix 
Hamiltonian Hf2_ is an odd function of k [see Eqs. (|21f) 
and H22|) ]. and (ii) the last two terms on the right-hand 
of (|25|l also disappear, for they describe the scattering in 
and out of a momentum state, hence the average of the 
total rate of change of occupation of a momentum state 
due to scattering is identically zero. 

To solve the kinetic equation for the electron spin den- 
sity S(i), we first need to evaluate the density matrix 
Pnki^) i-"^ ^Q- CZI)- In the first order approximation in 
A(n, k), the density matrix p^^it) is determined by the 
following equation: 

d 



dt 



fl2it) 



^nk ' Pnk V-l 



E 



dk' 

W) 



%w}!lo^y){^l{t)-iX,{t)), 



(29) 



which is obtained from Eq. ((SHJ- From Eqs. lj2II)-(|2nil it 
can be seen that the relaxation rate of the spin density 



-i \ttSO JO), . 



E 



dk' 

(M 



■w}^l{k,k')U;:l{t)-iiX{t) 



(30) 



E. QUASI-ELASTIC SCATTERING 



To make the problem more tractable and to express the 
solution to the spin density matrix equations [Eqs. I|27|l 
and (|30|l ] in a physically revealing way, the collision inte- 
gral on the right-hand side of Eq. (|30|) is further approx- 
imated by assuming the momentum scattering is quasi 
elastio^. In addition we assume that the energy spec- 
trum is isotropic (i.e., a spherical Fermi surface, eniy{k) 
— e„cr(|k|)), and this implies that |k'| = |k| as de- 
manded by energy conservation. Under these assump- 
tions, the momentum scattering rate Wrm' (k, k') appear- 
ing in Eq. I|20() depends only on the magnitude of k and 
the scattering angle ^kk' between k and k'. Thus the 
momentum scattering rate can be written as 



li:w(k,k') =!£„„' (k-k,fc). 



(31) 



where cos ^kk' = k • k . All spin coherence times will be 
calculated in this paper assuming quasi-elastic scattering. 
Using Eq. 1)31(1 in ((3()|l . the resulting simplified expression 
is given by 



—I 



dk' 



E 



r{5) 



(27r)3 = 



W:r^'(k-k,fc) 



p!^w 



il'it) 



(32) 



The coupled equations l|77|l and l|5^ are quite general 
and they can be applied to both three-dimensional (bulk 
semiconductors) and two-dimensional (quantum wells 
and superlattices) systems. We will treat these two cases 
separately in the subsequent sections. 



F. SPIN DECOHERENCE IN BULK 
ZINCBLENDE CRYSTALS 

In the case of bulk semiconductors, the coupled 
Eqs. (|27|l and H32|) are explicitly written as 



d_ 
dt 



m 



Tr 



E 



dkP d 



-T^^ 



E 



27r2 dt = 



dkk^ 



(Pnk(i))i 



27r2 



dOi. 



Att 






= E 



dk'k 



'l,'2 



dOc 



ku/('5) ft 



li:r^'(k-k,fc) 



27r2 J At: — "" 

X (ilit) - i'Ut)) ^ (34) 

where dO^ is the differential solid angle in k-space and 
(Pnk(O)k - {^i")"^ Jd0^^pnu{t). In order to solve 
the coupled equations (|S^ and l|M)l . we shall expand 
iifik ' ii^nn' (k • k , A;) , and p^j^ (i) in a series of spherical 
harmonics^ 



sf.£ = EM°("'^)^'™( 



Im 



^lit)^T.~^2(n,k,t)YU0,< 



(35) 



(36) 



^it!'(k • k , fc) = ^^it!.(fc)^/(cos^kk') 

= Elin'/W^^^r^i^'.-^Oii™!^,-^), (37) 



where Pi (^kk' ) is the Legendre polynomial and the second 
equality in Eq. (|37|) is obtained by using the addition 
theorem for spherical harmonics'*^. 

Substituting Eqs. lj23)-(|SZ|l into Eq. ^, Multiplying 
both sides of the resulting equation by Yi* ^, {9 , (f>) , and 
integrating over the solid angle dO^^ we have 

£lm("> k, t) = ^n{n, k) [1S0(„^ k), {gn^.{t))^ , (38) 

where the momentum scattering time matrix T;(n, k) is 
given by 

i!r'\n,k) = £„(fc)E/^E-'(k •£',&) 



l-PKk'-k) 



(39) 



In obtaining Eq. (|39|) . the delta function 
(5 (£„'o-'(k') — e„CT(k)) contained in the definition of 



the momentum scattering rate W£„„'(k, k') in Eq. l(^ 
is eliminated by the integration dk' . Here V (fc) — 

'Dna ik)I_, where Vna (k) is the three-dimensional 
density of states of the nth band with spin index a, and 
/ is the 2x2 identity matrix. Note that the density of 

states Vila {k) is the same for both up spin and down 
spin due to Kramers degeneracy: V^-. ' (k) — V^. (k). 



As the scattering potential is spin independent, the 
momentum scattering time matrix is diagonal with 
respect to the spin indices and can be rewritten as 
(33) I.i{n,k) = Ti{n,k)£, where Ti{n,k) is the momentum 
relaxation time. The I = 1 component of the momentum 
relaxation time Ti(n, fc) corresponds to the transport 
time Tti(n, fc) appearing in the mobility expression [see 
Eq. ^ in Secirrj]. 

Inserting Eq. H38|) into Eq. (|33|) and performing the 
integration over the solid angle dO^ and the sum over all 
states n, the resulting expression is 



l-S(f) ^Tr 

dt (27r)3?i2 



J2 f dkk^{~-l)'^Ti{n,k) 

_ Imn 



m„An.k), m-.nin,k),{p^^{t))^ 



where H_f^{n, k) is written as 



, (40) 



\ldO. 



f^(n,k)-oi 1,™(0,0).(41) 



To evaluate the trace on the right-hand side of Eq. I|40|l , 
we shall write the spin density matrix (/Onk(i))k in a gen- 
eral form as 



(P"k(i))k 



1 



g{n,k)l+h{n,k)S{t) a 



(42) 



where g(n, k) and h{n, k) are electron distribution func- 
tions which are determined by the initial condition at 
t = to [Eq. I|14|l ] and carrier statistics. In general these 
distribution functions can be written as^^^ 



liji.k) = 



f^+\n,k) + f^-\n,k) 



E/|^(/^^("'^) + /^'("'^)) 



, (43a) 



h{n, k) 



ft\n,k)-ft\n,k) 



E 



dkk^ 

27r2 



{f!,-^\n,k)-ft\n,k)) 



, (43b) 



where /q {n, k) are the Fermi distribution functions for 
spin-up and spin-down electrons with the corresponding 
Fermi energies E^ ' . 



After substituting Eq. l(12)) into Eq. l(^ and taking the 
trace on the right-hand side, we have the final form of the 
kinetic equation for the electron spin density, expressed in 
terms of the precession vector and momentum scattering 
time: 

^S(i) = -l^Y. J dke{-iy^T,{n,k)h{n,k) 

Imn 



[p.i-m{n, k) ■ rtimin, k)] S(<) 



rii^rnin, k) ■ S(i) ) ftim{n, k) 



where 



ni„,{n,k)= / d0^n{n,'k)Yi 



m\",VJ^ 



and 



fl(n,k) = ^n,™(n,fc)l^,„(0,, 



(44) 



(45) 



(46) 



Im 



In spherical coordinates, Eq. (|22|l can be written as 
Cl{n,k,9,(j)) — ~n{n,k,Tr — 9,tt + (p). Therefore, from 
Eqs. ^^ and ^^, we have 

Y,^im{n,k)Yi^{e,^) 

Im 

= y^ -Vtlra[n, k)Yira{-K - g, TT + 0) 

Im 

= Y,^hn(n,k)Yi^{e,(^){-l)'+\ (47) 



Ira 



where the last equality is obtained using y;m(7r — 0, tt- 
= Yi^{e,(j)){-1)''. Equation g7|) yields 

hi„,{n,k) = 0,,„(n,fc)(-l)'+i. 



(48) 



In order to satisfy this condition, Jl;.,„(n, fc) must vanish 
for even I. 

Rewriting Eq. (|44(l in tensorial form, we have 



d 



-5,(i) = -5]r,,5^(t), 



dt 



(49) 



where (i,j) G {a;,y, z} and the components of the spin 
relaxation rate tensor F are given by 



Yij = .g y^ y^ / dkk'^Ti{n,k)h{n,k) 

Imn 



X! ^*j'hnin, k)nj,im{n, k)Si. 



namin,k)n*i^{n,k) 



(50) 



In deriving the last expression, we used $!;_,„ (f^,fc) — 
{-l)"'n*^(n,k), noting that J7(n,k) is real. The lon- 
gitudinal Ti and transverse T2 spin relaxation times are 



now directly obtained from Eq. 150|) [see Sec. Ill L| for de- 
tailed discussion] . In our spin coherence time calculations 
for bulk semiconductors h{n, k) — fo{n, k) [1 — /o(n-, fc)] 
^j^(3D)-j-i^ where /o(n, k) is the Fermi occupation func- 
tion and n!^^^^ is the three-dimensional electron density. 



G. PERTURBATIVE EXPANSION FOR BULK 
ZINCBLENDE CRYSTALS 

In general the electron spin splitting must be calcu- 
lated numerically; however it was argued by D'yakonov 
and Perer«2Si22i2£ that, to a good approximation, the 
electron spin splitting for bulk zincblende crystals can 
be evaluated using a perturbative approach. In a third 
order perturbation approximation (the lowest order for 
non- vanishing spin splitting), the spin splitting is propor- 
tional to k^ and the components of the precession vector 

28.29.30 



^.(l,k) = 



f],(l,k) = 



^.(l,k) 



fi^ac 


{2mlE^)h 


h etc 


{2mlE,)i 


fl CXc 



T^xyky fc^l. 



ly^z 



-K 



{2mlE^)-^ 



T'^z\k^ ky], 



(51a) 



(51b) 



(51c) 



where rric is the electron effective mass. Eg is the band- 
gap energy, ac = jc^~^Eg' (2toc)^^^ is a dimensionless 
constant related to the spin splitting of the conduction 
band, and 7c is the conduction-band spin-splitting coef- 
ficient (see, for example, Ref. ,44, for calculations of 7c for 
zincblende semiconductors). This corresponds to taking 
only Z = 3 in the series expansion of Eq. H46|) . and the 
non-zero expansion coefficients are 



^y3±i{^,k) = inx3±i{l,k) = 

^y3±3{^,k) = ina:3±3{^,k) = 

o n i^^ 27rih^a,k^ 

"^3±2(l,fc) = 



iir^ti ark^ 



{42mlEg 



(70m3i?g)i 



{lObmlEg 



(52a) 
(52b) 
(52c) 



In this case, the kinetic equation [Eq. H50|l ] for the spin- 
polarization vector can be simplified to 

§-^S,{t) = -ljdkk^?,{l,k)h{l,k) 



E 



(fi,(l,k)17,(l,k))£5.(t) 



- (f],(l,k)l],(l,k))j^5,(i) 



(53) 



and the components of the spin relaxation rate tensor Ti 
in Eq. (|1^ are reduced to 

^JdkPh{l,k)?3{l,k) 

E,(^|(l,k))£-(f}?(l,k)}£l ii^ 



^Jdkk^h{l,k)?3{l,k) 



if* ^3, 
(54) 

For a Maxwell-Boltzmann distribution for h{l,k), 
D'yakonov and Perel'^SiSS obtained the foUowing analyt- 
ical expression for Ti: 

^ ^ Ea^TpikBT)^ 
Ti 



with 





n^Eg 


1 






16 
35 i 


,im^^^ 


> 


+ 


5 
2 



(55) 



(56) 



where ^b is the Boltzmann constant, Tp is the average 
transport time which is directly related to the electron 
mobility through /i = eTp/rric, and the numerical con- 
stants rp'^i and v depend on the momentum scatter- 
ing mechanism (and are given in Table P) . In the case 
of a degenerate electron density^^-'^^ , the expression for 
T\ above is modified by simply replacing k-^T with the 
Fermi energy E-p and setting S = 32/(105?7'^'^^^). This 
perturbative approach provides an adequate description 
for electron spin decoherence in bulk zincblende crys- 
tals. The comparison between the perturbative and non- 
perturbative approach is made in Sec. IIVI 



^fK,Pi?^w =E 



27r 



27r 



'k;, 



27r 



x^S'(K|| -Kji-i^lhi^.) [^I{t)-P%\t)) , (58) 

where K is the superlattice wavevector with K^ along the 
growth direction z of the superlattice, Kii is the in-plane 
wavevector, the quantum number L is the superlattice 
band index that enumerates all superlattice eigenstates, 
and (p£k(0)k = (^i")"^ /'^<^k denotes the averaging 
of PCK-it) over all directions of the in-plane wavevector 
Kii. The momentum scattering in this case is essentially 
two-dimensional, therefore the electron scattering transi- 
tion probability Wrj-iCK^^ ■ 'K'..,K\\,Kz) depends on the 
magnitude of the in-plane wavevector Ku and angle be- 
tween the initial and final in-plane momenta (Kii and 
K||) for each K^. Under this assumption, we are effec- 
tively solving the coupled Eqs. (|57|) and H58|) of an ideal 
two-dimensional system for each value of Kz and then 
averaging the final results over the momentum along the 
growth direction. This approximation is quite reasonable 
because the band structure along K^ of the layered struc- 
tures we will consider here is relatively flat (i.e., almost 
dispersionless)^2iiS. 

To solve the coupled spin density matrix Eqs. H57|l and 
•EHl' ^£K|| - ££k(0, and Kcl (K|| • K|| , X,, , K^) are ex- 
panded in a Fourier series in the Ky-plane for each value 
of K, ^2^ 



£^K„ = 2^ Sr (A ^11 , i^.) exp(z;0jt|| ). 
I 



(59) 



PcKit) = Y,-^,'\c,K\\,Kz,t)exp{U<l>^), (60) 



H. SPIN DECOHERENCE IN 
NANOSTRUCTURES 

In this section, we consider the electron spin dynam- 
ics for both quantum wells and superlattices simultane- 
ously. Quantum wells can be viewed as superlattices with 
large barrier thickness in which the wells are practically 
isolated; therefore, the spin relaxation in quantum wells 
and superlattices can be treated on the same footing. For 
two-dimensional systems, it is convenient to use a cylin- 
drical coordinate system. The coupled equations H27|l 
and (|32|l for the two-dimensional systems are explicitly 
written as 

dK, 

~2^ 



I^W 



Tr 



E 

£ 



dKi\K\i 



d 



27r 



-Tr 



E 



dKz 

27r 



^^(^£K(t))K„^ 

di^iiXii [d4>^ 



27r 



27r 



^£K;P£Ki*) ^ 



(57) 



^S'(KrKj|,i^l|,if,) 

I 

Substituting Eqs. ^^-^ into Eq. (gHll, Multiplying 
both sides of the resulting equation by exp(— j^t/ij^ ), and 

integrating over the dqit^ , , we have 
pp)(£,X||,if„i) 






lfO(/:,if||,if.),(££K(i))K„ 

(62) 

where the momentum scattering time matrix 
?/(£, K\\ , Kz) is given by 



li\C,K\\,Kz 



1 dS 



KnK', 



2-n 



UcY. 



x[i-cos(;Sk||K;|) 



Wcc(K||-K' if||,i^,) 



(63) 



where cosSkmK' 



K 



Kj|. Here 2c 



-n(2-D) r ■ 



the two-dimensional density of states of the Lth band 
with spin index a. Note that the constant density of 

states 2?}^^ is independent of the superlattice wavevec- 
tor and is the same for both up spin and down spin due 
to Kramers degeneracy: V)^^ ' — T)]^, ' . As in the three- 
dimensional case the momentum scattering time matrix 
is diagonal with respect to the spin indices and hence it 
can be written as z.i{C,K\\,Kz) = ti{C,K\\,Kz)L For 
superlattices the transport time rtr(/^, -ftTy, X^) from the 
two-dimensional mobility is given by jj^i{CtK\\,Kz) as 
in Eq. H63|l (see Sec. III.TJI . Notice the different angular 
dependence of the transport time in bulk crystals and 
superlattices. 

Substituting Eq. (|()^ into Eq. (|57jl and performing the 
integration over 



^k' 



we obtain 



l^<" 



-1 



(27r;i): 



rTr 






^^"(A if II, if.), ^^V(Aif||,if.),(P£K(t))K„ 



where Hf^iC, if||, Kz) is obtained from the relation 



(64) 






Kr{C,K\\,Kz) 



27r 



H 



so 



XK|| exp(,-i<(^j^ 



^j 2^ [f2(AK).^exp(-*;0K||)- (65) 

To evaluate the trace on the right-hand side of Eq. (|64|l . 
we write the spin density matrix {p cK.it)) ^ in the form 
of Eq. (gll) 



(P£K(t))K 



g{C, K\\ , Kz)L+hiC, K\\ , Kz)S{t) -^ 

(66) 

where g(£, K\\ , Kz) and ft,(£, K\\ , Kz) are the electron dis- 
tribution functions which are written similar to Eqs. (|43|1 . 
Substituting Eq. (|H()|l into Eq. f^^ and taking the trace 
with respect to the spin indices, we obtain the kinetic 
equation for the electron spin density for superlattices 



where the components of the spin relaxation rate F are 
given by 

r^^ - ji^J2 J dK. J dK\\K\\niC,K\\,Kz) 



xh{C,Kn,Kz) 



J2^}iiJ^,K\\,Kz)nriiC,K^\,Kz)5^, 



-n*iiC,Kn,Kz)nu{C,Ku,Kz) 



(69) 



The last expression is obtained using Cli-i{C, K\\, Kz) 
= Q*i{C, K\\, Kz). In our spin coherence time calcula- 
tions for quantum wells and superlattices h{C, Kn, Kz) 
= fo{C,K^i,Kz) [1- MC,K\\,Kz)] (n(3^))-i, where 
fo{C, K\\ , Kz) is the Fermi occupation function and n^^^^ 
is the three-dimensional electron density. 



I. PERTURBATIVE APPROACH FOR 
QUANTUM WELLS 

The extension of the perturbation approximation in 
bulk crystals to two-dimensional systems is originally due 
to D'yakonov and Kachorovskii^i. In a third order ap- 
proximation, the precession vector for structures grown 
along (001) direction is^ 

n(l,K||) 

= nxi{i,Ki\)cos{(j)^^^) + nx3{l,K\\)cos{3(j)^^^), 

hy,{l, if||) sin(0ji|| ) + f^,3(l, if||) sin(30K„ ), 

(70) 

with the angular components Qji{l, Ki^\) given by 



^.i(l,if|| 



fe2 

-Ti Uc 



-r!^i(i,X|i 

IrUcEia 



8{2mlE^ 



Ki\K 



h 



(71a) 



17,3(l,if||) = fi,3(l,if||) 



-h^ 



3.^^^^ (^ib) 



8(2™3£;g)2 



n^i{C,K,uKz)-ni{C,K«,Kz)]S{t) 



{n_i{C, K\\,Kz)- S(t)) ni{C, K\\,Kz 



(67) 



Rewriting Eq. 1)67(1 in a matrix form equivalent to 
Eq. l(^ . we have 



d_ 
di 



s{t) = -rs(t), 



(68) 



where E'lc is the electron confinement energy of the first 
quantum-well state. In this approach, the crystal effec- 
tive magnetic field along the growth direction vanishes 
because the z-component of the precession vector is iden- 
tically zero. Using Eqs. H71|l in Eq. (|69|l . the components 
of the spin relaxation rate tensor are 

^ XX ^ ^ yy ^ T^i zz 
3 

= iE /rfif||if||/i(l,if||)?Kl,if||)^'i(l,if||)- 
1=1 •' 

(72) 
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In Ref. Is^, the Ku term in Eqs. ifTT)) was considered 
to be small compared to the linear K\\ term, and conse- 
quently it was dropped, as well as the I = 3 term. Under 
this assumption, D'yakonov and Kachorovski:^^ obtained 
an analytical expression for the Ti of a non-degenerate 
electron gas after averaging Eq. H72|l over the Maxwell- 
Boltzmann distribution 



1 



2alElTpkBT 



h^E, 



(73) 



and 



y'rfi^„4Mi,i^ll)n(i,i^ll) 



dKi\Kf^h{l,K\\ 



(74) 



where Tp (see, for example, Ref.|4J|) is the average trans- 
port time, related to the electron mobility via /i = 
eTp/rric [see Eq. (IHTll in Sec. Ill ,Tj . In the case of a de- 
generate electron gas, Eq. (|73|l is modified by replacing 
k^T with the Fermi energy Ep. The validity of these 
approximations is examined in Sec. IIVI 



J. TRANSPORT TIME AND MOBILITY 

To evaluate the momentum relaxation time Ti{n,k) 
for bulk semiconductors [see Eqs. (|39|l ]. we first extract 
the transport time Ttr{n^k) [i.e., Ti{n,k)] from the ex- 
perimental measured mobility by solving the Boltzmann 
equation using a fourteen-band non-parabolic electronic 
structure for the bulk semiconductor (see Sec. lIIip . Once 
the transport time is determined, the momentum relax- 
ation time Ti{n, k) for / > 1 is expressed in terms of the 
transport time Ti(n, k). The Boltzmann equation is writ- 
ten in the standard form (see, for example, Ref. l35j) 



dfjr, 71,^,1) 
dt 



v(n,k)- Vr/(r,?i,k,i) 



eh-^Fir,t)- Vk/(r,n,k,t) 



,k.t 



[/](75) 



where e is the elementary charge, v(n, k) is the electron 
velocity in the nth band, F(r,i) is the electric field, and 
the collision integral on the right-hand side of Eq. (|75|l is 
given by 



/r,„.k..[/] = -E/7|^^""'(k',k) 
X [/(r,n,k,i)-/(r,n',k',t) 



(76) 



In the relaxation time approximation, the collision in- 
tegral takes the form 



Ir,nM,t [/] 



/(r,n,k,^) -/o(r,ri.,k,t) 

Ttr("-, k) 



(77) 



where /o(r, n, k, t) is the thermal equilibrium distribution 
function and rtr(n.,k) is the transport time. In steady 



state, for a spatially homogeneous and time-independent 
electric field, the solution of the Boltzmann equation H75|) 
linearized with respect to the electric field F is written 
as 

/(n,k) = /o(n,k)-e?i-Vtr(n,k)F. [Vk/(n,k)]. (78) 

The electron mobility tensor /x is defined as 

Vdrift = M • F, (79) 

where Vdrift is the average drift velocity and it is obtained 
by averaging the momentum-dependent electron velocity 
over the non-equilibrium distribution function 



Vdrift 



Y^ f dkf{n,k)v{n,k) 

n •' 

J2 fdkfin,k) 



(80) 



Thus the electron mobility tensor for bulk semiconduc- 
tors can be obtained from Eqs. (jZHJ, IZHI), and ((HOI, and 
the final expression is written as^^*^ 



/^ 



Bulk 



l(30)?i' 



■E 



dk 



-Ttr{n,k) 



dfa{n,k) 



Vk-B„,(k)|[Vk^„<x(k) 



dE^^ik) 
(81) 



where n^^^> is the three-dimensional carrier density and 
Ena{k) is the electron energy in the nth band. Equa- 
tion (|81|l connects the transport time rtr(n, k) to the 
experimental measured mobility /x^""'. For zincblende 
semiconductors, the mobility tensor [Ea.f I81|l ] is diago- 
nal with identical elements owing to the cubic symmetry 
of the zincblende materials. 

The electron mobility tensor yF"^ for superlattices is 
obtained in a similar fashion, solving the linearized Boltz- 
mann equation within the framework of a generalized 
K-p envelope function approach (see Sec. IIII|I . In this 
case, the electron mobility tensor for superlattices takes 
the same form as that for the bulk case, Eq. (|81|l (see 
Sec. nil El for a brief discussion on symmetry). For 
zincblende type superlattices, the cubic symmetry of the 
constituent semiconductors is preserved in the in-plane 
direction, but it is broken along the growth direction. 
As a result, the diagonal element of the mobility tensor 
along the growth direction takes a different value than 
those along the in-plane directions. 



K. ORBITAL SCATTERING TIME 

To calculate the orbital relaxation time [see Eq. (|S^ ]. 
the functional form of the orbital scattering rate 
Wnn'ik 'k! ^k) is taken from standard expressions for 
ionized impurity (II), neutral impurity (NI), or optical 
phonon (OP) scatteringSSiffi. The scattering time Ti{n, k) 
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for Z > 1 is expressed in terms of the transport time 
Ttr(n, A) — Ti(n, fc), and the ratios for ionized impurity, 
neutral impurity, and optical phonon scattering are sum- 
marized in Tabled The transport time Ttr(n, fc) is ob- 
tained by inverting the Boltzmann transport equation 
()81|l using a fourteen-band nonparabolic electronic struc- 
ture for the semiconductor (see Sec. lIII|l . In Eq. H81|) . the 
transport time is expressed as a function of energy and 
it has the form of rtr(n, k) oc E^^(k). The values of i^ for 
various scattering mechanisms are listed in Table Q] 

The approach to evaluate the momentum relaxation 
time Ti(C, Ku, K^) for two-dimensional structures is very 
similar to that of the bulk case. The functional form 
of the scattering rate W^££'(Kj| • Kj,, ^n, X^) is again 
taken from standard expression for ionized impurity, 
neutral impurity (such as arises from superlattice in- 
terface roughness), or optical phonon scattering. The 
Ti{C, K^^jKzYs differ for different mechanisms and the 
relationship between the momentum relaxation time and 
the transport time for various scattering mechanisms 
is tabulated in Table ^1 The connection between the 
momentum relaxation time and the transport time in 
two-dimensional systems is very simple, and the ratio 
of Ti{C, K^^,Kz) and Tti-{C,K^^,Kz) can be expressed as 
a function of the angular index I of the field compo- 
nent. The effective time for field reversal ti{C,K^^,Kz) 
depends on the angular indices I of the field compo- 
nent Qji{C, K^^, Kz). For example, an Z = 1 compo- 
nent il,ji{C,K^^,Kz) requires a 180° change in the an- 
gle of K to change sign, whereas a.n I — 3 compo- 
nent Qj3{C,K\\,Kz) only requires a 60° change, so typ- 
ically T3{C,K\^, Kz) < Ti{C, K\^,Kz). The momentum- 
dependent transport time has the form of Tti-{C, K\\ , K^) 
oc £'£g(K), and the values of v for ionized impurity, neu- 
tral impurity, and optical phonon scattering are given in 
Table ini 



rp-l 



T^E 



dk 



-i7i(n,k)+i/|f(n,k) 



X r(n, k)/i(n, k). 



(82b) 



where the electron distribution function h{n, k) is defined 
similar to Eq. I|43b(l . T = (g/ZB^~ )^, A*b is the Bohr mag- 
neton and g is the electron g-factor. In a crystal with in- 
version asymmetry and spin-orbit coupling there is a spin 
splitting described by the Hamiltonian given in Eq. H21|l . 
The momentum-dependent effective magnetic field due to 
inversion asymmetry is H(n, k) — T^^/^r2(n, k). As the 
electron is scattered from k to k' via ordinary orbital (not 
spin-dependent) elastic scattering, the effective magnetic 
field H(n, k) changes direction with time. If the crys- 
tal is cubic, then i?^(n, k) = Hy{n,k) = iJ|(n, k), so 
Hl{n,k) = 27J|(n,k) and T2 = Ti (see Ref.0). 

The relationship between Ti and T2 differs, however, 
for systems of lower symmetry, such as quantum wells or 
superlattices. For a (OOl)-grown quantum well, one finds 
in perturbative expressiongSSii^ that the fluctuating field 
along the growth direction vanishes. In this case^ 



rl 



T^\a) = Tj;\a = 0) -(1 + cos^ a) 



T-\a)=T-\a = Q) 



-(2 



sin a] 



(83a) 



(83b) 



where a is the direction between Ho and the growth di- 
rection. Thus T2 ranges from 2Ti/3 to 2Ti depending on 
a. By similar arguments, for (llO)-grown quantum wells 
the effective crystal magnetic field is entirely along the 
growth direction, and 



T-\a) =Tfi(a = 0)si 



sm a. 



(84a) 



L. GENERAL FEATURES OF THE 

LONGITUDINAL Ti AND TRANSVERSE T2 SPIN 

RELAXATION TIMES 



This section considers both bulk and quantum- well sys- 
tems. In the motional narrowing regime the electronic 
spin system is subject to an effective time- dependent, ran- 
domly oriented magnetic field H(n, k) which changes di- 
rection with a time rin, k) that is much shorter than the 
precession time of either the constant applied field Ho or 
the random field. The coherence times in bulk zincblende 
semiconductors depend on the transverse 7Jj^(n, k) and 
longitudinal Hu {n, k) components of the random field, 
according to^ 



Tfi 



T 






dkHl{n,k)T{n,]<.)h{n,k), (82a) 



T2-1 (a) = Tf 1 (a = 0)[i(l-fcos2 a) 



(84b) 



Thus although T-^ ^{a = 0) vanishes, the same is not the 
case for T~'^{a) for any a. From Eqs. ^, ^, ifS^ . 
and H84|l . one can immediately identify 



rfi(o) = r. 



T, 



-\o) = i(r.. + r,,). 



(85a) 



(85b) 



M. CONNECTION BETWEEN Ti AND Ts TIMES 

In this section, we briefly discuss the relationship be- 
tween the longitudinal spin relaxation time Ti and the 
spin-flip time Ts, which is usually refered to as the "spin 
relaxation time" in optical experimentsii*i2ii^iiSiiS*ii. 
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In spin-dependent pump-probe experiments, a non- 
equilibrium spin population is created by circularly po- 
larized light, and the spin orientation of the excited elec- 
trons depends on the helicity of the incident light. Once 
the unbalanced electron spin population is excited in the 
conduction band, it relaxes toward equilibrium both via 
carrier recombination and via flipping spins. It is usu- 
ally assumed that the unbalanced spin population is de- 
scribed by the rate equationsii^ 

dN+{t) _ N+it) N+{t) N^{t) 

— -- — — 1 , (soa) 



dN_{t) _ N^{t) N+{t) N^{t) 



(86b) 



where N^(^^){t) is the number of circularly polarized elec- 
trons in spin-up (down) state, r^ is the electron-hole re- 
combination time, and r^ is the spin-flip time22ii&. Tak- 
ing the difference between Eqs. H86a() and l)86b|l gives 



dD{t) 
dt 



1 2 , , , 



(87) 



where D{t) is defined as D{t) = (iV+(t) - iV_(i)), and is 
directly proportional to the macroscopic magnetization 
of the system. The solution of Eq. H87() takes a simple 
exponential form 



A. FOURTEEN-BAND RESTRICTED BASIS 

The Hamiltonian H pertinent for the bulk constituents 
of heterostrutures is given by Eq. Q. It is solved within 
a fourteen-band restricted basis set of the constituents, 
consisting of two conduction antibonding s states (s, 
or the Eg representation), six valence bonding p states 
(Eg © E7), and six antibonding p states (p, also Eg © E7) 
(see Table IIIip ^^i^^i'^''' . These fourteen basis states are 
characterized by the total angular momentum quantum 
number (J, J^). The total angular momentum is J = L 
+ S, where L is the orbital angular momentum and S 
is the spin angular momentum. The orbital angular mo- 
mentum takes the values of L = for the s and s states 
and L = 1 for the p and p states. This is the minimum 
number of basis states that are required to generate spin 
splitting non-perturbatively (see, for example, Ref. 41). 

The electronic band structure of bulk GaAs at room 
temperature is shown in Eig. Q^a) together with its pa- 
rameters defined in the k-p formalism. Here E^ is 
the fundamental energy gap [between ^(Eg) and ^(Es) 
states], Ec is the energy difference between ^(Ey) and 
p(Eg), A is the spin-orbit splitting of the p states, and 
Ac is the spin-orbit splitting of the p states. The three 
interband momentum matrix elements (Pqj P\^ a-nd Q) 
that couple the states away from the zone center in the 
k-p model are 



iPo = (.s\u^-) - (5|p,|y^') = (^Ip.iz''), 



(89a) 



D{t) = L>o exp {-tlTr) exp {-^tJT,), 



where Dq is the initial excited electron spin population. 
Eor an electron-hole recombination time much greater 
than the electron spin relaxation time (t^ S> t^), the 
unbalanced spin population decays exponentially toward 
equilibrium with a decay constant Ts/2. Thus the macro- 
scopic magnetization of the system relaxes exponentially 
with a time constant Ts/2, and the longitudinal spin re- 
laxation time Ti = Ts/2 (see, for example, Refs. |33 and 

SF 



III. FOURTEEN-BAND K p THEORY 

Calculations of the longitudinal time T\ and the trans- 
verse time T2 require knowledge of the precession vector 
J7(£, K), which is directly related to the energy eigen- 
states and eigenvectors of the system. In this section, 
a generalized fourteen-band non-perturbative theory is 
presented for the electronic band structure of zincblende 
semiconductor heterostructures based on the superlat- 
tice K-p theory22i22iSli2SiS&, which is a generalization of 
the bulk k-p theory for systems with periodicity on the 
scales larger than the bulk lattice constant. This peri- 
odicity exists only in the superlattice growth direction, 
which in this case is taken to be the z direction. 



zPi = (5|p,|x^) = (5|p,|y^) = (^Ip.Iz^), 



(89b) 






{Z"\Py\X^) 



{Y^lp.lX'^) 

{x^\p.\yi- 



The energy splitting at finite crystal momentum k of 
GaAs shown in Fig. ^a) is a result of spin-orbit inter- 
action and inversion asymmetry in zincblende semicon- 
ductors. At zone center, the eigenstates are mixtures 
of different spin eigenstates but are still doubly degen- 
erate, and so a "pseudospin" degeneracy can be identi- 
fied. Away from zone center, the energy eigenstates of the 
crystal Hamiltonian H are certainly not spin eigenstates, 
and it is not even possible to choose a pseudospin ba- 
sis which diagonalizes the Hamiltonian for all directions 
of k. Tables ]V\ and IVll present the parameters for the 
systems described in Sec. IIVI which are obtained from 
fits to experimental measurements of the conduction- 
band mass, the heavy- and light-hole masses, and the 
g-factor given in Refs. lisL l49l and |53- The parame- 
ters in Tables IVl and IVll obtained from these experimen- 
tal observables have been used in several recent publica- 
tions describing various optical and electronic properties 
of semiconductors^ii^SiS, and have not been altered for 
our spin coherence calculations in bulk and quantum-well 
semiconductors. 
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B. SUPERLATTICE WAVE EQUATION 

The Schrodinger equation for heterostructure superlat- 
tices is written as 

H'''^{r\C,S,K) ^ Ecs{K){r\C,S,K), (90) 

and the superlattice Hamihonian H^^ is written as 

-'Via 



i?SL = J2 HA{v). 



j=l 



.(r) = 



1 if r G ith layer, 
if r ^ ith layer, 



(91a) 



(91b) 



where Hi is the crystal Hamiltonian of the ith layer 
[Eq. ||2Jl], Mayor is the number of layers in the super- 
lattice unit cell, |£, 5, K) is the superlattice eigenstate 
for a carrier with wavevector K and spin S € {T? i} in 
the £th band, and Ecs{K.) is the corresponding super- 
lattice eigenenergy. To obtain \C,S, K) and Ecs{^), we 
first solve the Schrodinger equation [Eq. (|90|l ] at the zone 
center in which the superlattice wavevector K is identi- 
cally zero. 



ZONE-CENTER SUPERLATTICE STATES 



A zone-center superlattice wavefunction is expressed in 



terms of bulk electronic structure parameter: 



{r\N, S,0) = }^ FNSnair){r\n, a, 0), 



.22.26 



(92) 



where \n, a, 0), defined in Table lTTll are the corresponding 
zone-center Bloch states of the constituent bulk semicon- 
ductors, the expansion coefficients F]ysna{T^) are slowly- 
varying envelope functions on the scale of the constituent 
bulk semiconductor lattice constant, the quantum num- 
ber N is the superlattice zone-center band index, and S is 
the spin quantum number that enumerates both the spin- 
up and spin-down states. The ordering of the basis states 
given in Table Hill is such that the matrix Hamiltonian is 
block diagonal at k = 0. Here, for convenience we adopt 
the convention that lower and upper cases label bulk and 
superlattice states respectively. Equation (|92|1 implies 
that the zone-center superlattice states (r| A^, S, 0) are not 
spin eigenstates due to the mixing of different angular 
momenta. Therefore, the zone-center states (r|A^, |,0) 
and (r| A^, j, 0) are generally referred to as pseudospin-up 
and pseudospin-down states, respectively. By substitut- 
ing Eq. Ip?^ into l|?^ , Multiplying from the left-hand side 



of the resulting equation by the corresponding wavefunc- 
tion (n, cr, 0|r) of the zone-center Bloch state and inte- 
grating over a unit cell of the constituent bulk semicon- 
ductor, we obtain the zone-center Schrodinger equation 



J2 H^L'.'{r)FNSn'.'{r) = ENs{0)FNSna{r) (93) 



with 



nan' a' \ I 



2me 



'^'nn'^aa' 



£;^;'^,,(r)-f— p.Pr„^,-(r),(94) 



)Bulk 



^,^,(r) = (n,CT,0|p|n',f7',0), 



and 



£ 



Bulk 



(r) = (n,cr,0|iJ|n',cr',0) 



(95) 



(96) 



Note that within the total angular momentum basis 
states (i.e., zone-center eigenstates given in Table IhH . 
the bulk energy term _E^""^(r) in H^^{y) is a diagonal 
matrix whose diagonal terms are the eigenenergies of the 
constituent bulk states at zone center and they are 

E^f:t'A^)^E^f\v)5^^,5,,,, (97) 

where Ef^^^^{r) = £;,^""'(r). 

The presence of strain in superlattices due to lattice 
mismatch modifies the electronic band structure and 
can have dramatic effects on the electron spin coher- 
ence time. Consequently, it is important to take into 
account the strain induced change in the electronic band 
structure of superlattices. The implementation of the 
lattice-mismatchcd-induced strain within the framework 
of generalized K-p formalism is accomplished by mini- 
mizing the strain energy in the superlattice under the 
condition that the interfaces of the system are pseudo- 
morphia^. Consider a superlattice with periodicity Dsl 
along the z-axis. Projecting the full Hamiltonian ^ onto 
this fourteen-band restricted basis set with strain effects 
due to lattice mismatch taken into account, the effective 
strained k-p matrix Hamiltonian of the bulk constituents 
for (OOl)-growth superlattices takes the form o^ 



H 



Bulk 



ji I Fl ^ 
H^ 'h. 



(98) 



The sub-matrices I£in) and H^ are given respectively 
by ~ ~ 
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^T(i) = 



t^Bulk 

^cT(i) 



Pov, 



OVz 





PlVz 



-i\^PiVz 



cBulk 






hQvz 



^Povz 



cBulk 

^'T(i) 












cBulk 

'^'*T(i) 



QVz 



cBulk 

^i(T) 






IPlVz 
^QVz 








cBulk 

'T(i) 



kPiiz 

'^Vz 






cBulk 

^sT(i) 



(99a) 



H+ = 



—I 

i 

—I 

—i 

i 
—I 





IPov, 
\P,v- 

\PlV- 
\PlV- 



hPnv+ 








\Qv- 



kPoi- 






hQv+ 







3 





qQV- 

\Qv+ 





IPlV- 



\Qv- 












(99b) 



rBulk 



where Vz = rn^ ^pz, v± = m^ ^{px ± ipy), and £^ 
E^^^^ + ^,^t.ain + (2me)-ip^ and E^l""'" describes the 
energy shifts due to strain effects. The strained energy 
shifts are given bj*^ 



TT^Strain c 

■^cT(i) - Y' 



E 



Strain 



E 



Strain 



2' 



A 6^ 1 



(100a) 



(100b) 



S2Ni 



nil) --y + f +2(A^ + A^. + -A5^)5, (100c) 



E: 



' m)" = t + T " ^(^' + ^^" + i^^')^' (lood) 



with 



6c = 2Vai2ei\+e^), 



S, - 2Vb{e« - e^), 



an 
£11 = ^-1, 



(101a) 
(101b) 
(101c) 



e± 



a± 



1 



) C12 



- 1 



(lOld) 



(lOle) 



where e|| and e± are components of the strain tensor de- 
scribing the in-plane and growth-axis strains of the super- 
lattice, respectively, Oo, an, and a± are the unstrained, 
in-plane strained, and growth-axis strained lattice con- 
stants of the constituent semiconductors, respectively, 
cii and Ci2 are the elastic constants, and Va and VJ, are 
the hydrostatic and biaxial deformation potentials, re- 
spectively. The in-plane strained lattice constant an = 
Q'sub where Ogub is the lattice constant of the substrate. 
The energy shifts of the conduction bands are due to 
hydrostatic strain, while the energy shifts of the valence 
bands are due to biaxial compression [c.f., Eq. HlOlafl and 
Eq. IjlOlbjl ]. In the present treatment, the strain effects 
on the p states and the momentum matrix elements Pq, 
Pi , and Q are neglected. 

Note that f£| = I£i, i.e., the opposite-spin eigenstates 
of the constituent bulk states at zone center are two-fold 
degenerate, and ^_ = M}+- From Eqs. (EHl and l|^ 
it is easy to see that the zone-center matrix Hamilto- 
nian _ff ^*-^ (r) is block diagonal; therefore, the Bloch states 
\n, t,"0) and \n, J,, 0) are not mixed at zone center (see Ta- 
ble ^Oj. The zone-center eigenvalue problem [Eq. (jHSU 
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is solved in Fourier space in a similar method to that of 
Winkler and Rossler— . The spatially dependent terms 
in Eq. I|93|l can be expanded in terms of Fourier series 



^26.27 



H, 



zc 



/a'{z)= XI ^w<tV exp(i^gz) 



(102) 



FNSna{z)^ X FNSncrtxexp{ifiqz), (103) 

/^— — oo 

where q = 2tt/Dsl- By substituting Eqs. (|102|1 and Ij 103(1 
into ((93|l . Multiplying from the left-hand side of the re- 
sulting equation by exp(— i27ri^(72;), and integrating over 
the superlattice period, we obtain the zone-center super- 
lattice eigenvalue equation in Fourier space as 

where H'^^^,^, is given by the Fourier series expansion 



rVZC 



dzH^an'a'{-A^M~^m-A- (105) 



Diagonalizing the matrix H^^ in Eq. (|104|l yields the 
zone-center eigenvalues and the corresponding eigenvec- 
tors. Once the eigenvalues and eigenvectors are found, 
it is desirable to characterize the zone-center states in 
terms of the fourteen bulk characters designated by C|, 
hi, h, so|,7ix,7|, so|, C|, /i|, li, sox,7i|, 7^, and so^ (see 
Table Cnil. Noting that 



dr{N,S,0\r){r\N,S,0) 



VSL 



= V/ drF;js^^{r)FNSnair) 



'Vs 

x(n,cr,0|r)(r|n,CT,0}, (106) 

the fraction of each bulk character of the zone-center su- 
perlattice states is given by 



drF, 



NSna 



{r)FNSncrir) 



VsL 



E 



F 



NSnafi 



(107) 



where Vsl is the volume of the superlattice unit cell. The 
equality in Eq. (|1U7|I is obtained by using Parseval's the- 
orem for the discrete Fourier transform to replace the in- 
tegral on the left-hand side by the sum over the norm of 
the Fourier cofhcicnts on the right-hand side. Therefore 
the superlattice states at zone center can be classified 
according to the bulk character of their largest vector 
component, and the character Xn of the zone-center su- 
perlattice states is given by 



7C 

Xns = max 



/ , \FNSna-ti 



L M 



(108) 



where max[^„cr] stands for the maximum of the set 



D. SUPERLATTICE STATES AT FINITE K 

Once the zone-center superlattice eigenstates and 
eigenenergies have been determined, the superlattice 
states at K ^ can be obtained by application of su- 
perlattice K • p theorjiS^. A superlattice state \C,S,'K) 
at finite K is expressed in terms of the zone-center su- 
perlattice states 

(r|A5,K)=exp(zK-r)XC£5Jvs(K)(r|7V,5,0), 

NS 

(109) 

where CcsNsi^) are the expansion coefhcients and S £ 
{t, i} is the spin indices. Note that the superlattice in- 
plane wavevector Kii = (Kx, Ky) is equivalent to the bulk 
wavevector k|| = {kx,ky) in the in-plane direction due 
to in-plane translational invariance in one-dimensional 
superlattices, i.e., k|| is a good quantum number. 

By inserting Eq. (|109|) into (jHOI, Multiplying 
the resulting equation on the left-hand side by 
exp(— K • r)(A'', S', 0|r) and integrating over a superlat- 
tice unit cell, we have 



(110) 



N'S 



with 



^NSN'S'i^) 



2mp 



ENSN'S'iO) SnN'SsS' 



+ K-PjvSAr'S'(O), 

rup 



(111) 



where the zone-center matrix elements of momentum and 
energy are respectively given by 

Pnsn'S'{0) = (iV,5,0|p|iV',5',0), (112) 

Ensn'S'{0) = {N,S,0\H\N',S',0). (113) 

Making use of Eq. ((92|l . the momentum matrix elements 
connecting the zone-center superlattice states can be ex- 
pressed as 

{N,S,0\p\N',S',0) 



1 



— F 



3L - , J Vq, 

naiv a' ^^ 



dv 



no If o ■ 

FNSn'a'{Y)\n\a',Q) 



{n,a,Q\F*j^s^„{y) 



1 



— y 



dr 



Vsl 



F*NSnA^)FNSn'a'{Y) 



x(n, CT, 0|p|n',fT',0) 



-— y 



dv 



F, 



NSna- 



{r)pFi 



NSn'a' 



(r) 



x(n, CT, 0|n',cr',0) , 



(114) 
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The first term in Eq. H114() is evaluated in real space 
as its integrand is the product of three spatially depen- 
dent functions. The evaluation of the second term is per- 
formed in Fourier space to circumvent the problem of en- 
velope function variation on scales much smaller than the 
bulk lattice constant. Using the relaxation in Eq. (|l()3|l . 
the second term is simplified to 



1 



— y 



dr 



VSL 



F. 



NSjuj 



(r)pF, 



NSn'cr 



(r) 



{n,a,0\n',c7',0} 



M }, fJ-FMSnan^NSnati- 



nafi 



(115) 



Note that the first term in Eq. H114() is larger than the 
second term by a factor of TOe/mc, where rric is the zone- 
center conduction-band masa^. 



E. EVALUATION OF SPIN PRECESSION 
VECTOR 



e 



cos^-) =J2cl^^^iK)Cc^N^iK), 



N 



Z2 ^CiNI; (K)C£XAr| (K) 



(120b) 



(120c) 



N 



From Eq. (|120() , the trigonometric functions of the spher- 
ical angles (8, $) for the a;j/z-components of fl{C, K) are 
respectively written as 



sin(9) cos(<i>) = Re 



2E^^i^T(K)C£TA'T(K) 

(121a) 



N 



sin(9) sin(<i>) = Im 



2 22 ^CiN]' (K)C£|Ar| (K) 

. N 

(121b) 



To calculate the spin precession vector, it is convenient 
to introduce the pseudospin state |n;£,5, K) at finite 
momentum K with spin projection along the direction n 
parallel to the growth direction z (i.e., with quantization 
axis along the growth direction), 

\n-C,S,K) ^Y.^^SNs(K)Sss\N,S,0). (116) 

NS 

Any superlattice energy eigenstates |£,<S, K) can be ex- 
pressed as a linear superposition of the pseudospin-up 
|n;>C, t,K) and pseudospin-down |n;>C, J,,K) states and 
vice versa 



\h;C,S,K) = AslClK) + Bs\C,l,K) (117) 



with 



A^ = Bi = cos I — 1 exp I — — 



IT =^i 



Bi 



-A* 



sm 



e 



exp 



i(f> 



(118a) 



(118b) 



where and $ are respectively the polar and azimuthal 
angles that characterize the unit vector H(£, K) of the 
momentum-dependent effective magnetic field in the xyz- 
coordinate system. The electron spin precession fre- 
quency (Larmor frequency) is proportional to the spin 
splitting of the superlattice eigenstates and it is given by 



cos(e) = Y. [ I^^t^t(K)|' - [Ccimm\^ 

N 



(121c) 



where Rc[A] and Im[A], respectively, stands for the real 
and the imaginary parts of A. 

Taking into account the cubic symmetry of the 
zincblende constituent semiconductors, the evaluation of 
the precession vector n(£,K) for (OOl)-growth super- 
lattices can be simplified. The symmetry elements per- 
taining to (OOl)-oriented structures are: One fourfold 
rotoinversion axis (i4r2|j[001]), two twofold rotation axes 
(/'^^.[[[lOO] and ii'y||[010]), and two reflection planes con- 
taining the [001] axis and bisecting the angles between 
the [100] and [010] axes^^. All eigenstates with wavevec- 
tor K directed along the growth direction are twofold 
degenerate, and hence there are no spin splittings along 
this direction. Therefore, the evaluation of the in-plane 
precession vector r2||(£, K) can be simplified. One needs 
only evaluate ri^(£,K) and ^ylC.'K) over one eighth 
of the Kx-Ky plane (e.g., in the domain $ = [0, 7r/4]), 
and the rest can be obtained by symmetry consideration. 
This is accomplished by breaking the in-plane precession 
vector into eight quadrants on the Kx-Ky plane and by 
mapping each quadrant back into the first quadrant using 
the symmetry properties of the structure. The symmetry 
operations of the in-plane vector n||(£, K) are summa- 
rized in Table lIVl 



n{C,K)^h-^ Ec^{K)-Eci{K] 



(119) 



From Eqs. H116|) and (|117|l . we obtain the following rela- 
tions: 

\eMm^HQ)^Y.^hm{^)Cc^m{Vi). (I20a) 



N 



IV. RESULTS AND DISCUSSION 

In this section the electron spin decoherence in bulk, 
quantum-well, and superlattice zincblende semiconduc- 
tors is studied by examining the electron spin coher- 
ence times Ti and T2 using the fourteen-band non- 
perturbative theory developed in the previous sections. 
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These theoretical calculations are compared with avail- 
able experimental results. The fourteen-band k-p param- 
eters of the material systems used in these calculations 
are summarized in Tables IVl I VII and IVII 



A. SPIN COHERENCE TIMES IN BULK 
ZINCBLENDE SEMICONDUCTORS 

The understanding of electron spin dynamics in bulk 
semiconductors is essential because they are the basic 
building blocks of spin-dependent devices. To explore the 
electron spin properties of bulk zincblende semiconduc- 
tors, we calculate the electron spin coherence times using 
Eqs. l|5^ and (j82bp derived in Sec. Ill LI All the calcula- 
tions of electron spin coherence time for bulk semiconduc- 
tors are performed assuming ionized impurity scattering 
unless otherwise noted. 



1 . Ill- V SEMICOND UCTORS 

The electron spin coherence times T2S as a function 
of temperature T for GaAs are shown in Fig. Ola) , and 
the measured mobilities^ are displayed in Fig. OJb). 
Squares represent the experimental results by Kikkawa 
and Awschalom^. The agreement with experimental 
measurements for GaAs at the higher temperatures is 
goodSS, whereas for low temperatures other spin relax- 
ation mechanisms, namely the Elliot- Yafet (EY) mech- 
anism, may dominatc^^. At room temperature, the cal- 
culated spin coherence time T2 is approximately 216 ps 
with an electron mobility of 5.4 x 10^ cm^/Vs. From 
Fig.[2Ia) and|2Ib) the electron spin coherence times (both 
Ti and T2) are inversely proportional to the electron mo- 
bility [see Eqs. (|82a|) and (|82b|l ] . and hence the electron 
spin coherence times are inversely proportional to the 
momentum scattering times of the electrons. This im- 
plies that cleaner samples have shorter electron spin co- 
herence times. For semiconductors with cubic symmetry, 
Ti = T2 (see Seclim. 

Shown in Fig.|2c) are r2's as a function of temperature 
for GaSb and the measured mobilities'^ are shown in 
Fig. |2Id) . The shorter spin coherence times for GaSb are 
due partly to the larger conduction-band spin splitting 
Aii', which originates from a larger ratio R of the spin- 
orbit coupling A to the band gap Eg for GaSb (see Ref.l41l 
for perturbative expansions of spin splittings). The band 
structures of GaAs and GaSb at 300 K are shown in 
Figs.|Il;a)-(b), and i?GaAs = 0.24 and EcaSb = 1.05. 

For bulk zincblende crystals, all the electronic states at 
a given k are doubly degenerate along [100] directions [see 
Figs, n^a)-(b) for k||[100]], and the spin splitting along 
[111] directions is usually very small (not shown) while 
the spin splitting in [110] directions (see Figs, ma)-(b) 
for k||[110]) is the largest^. Figure|3Ia) shows the lowest 
conduction-band spin splitting as a function of crystal 
momentum (k||[110]) at 300 K, and Fig. |3fb) displays 



the electron spin splitting as a function of excess energy 
for k in the same direction as in Fig. Ola). For small 
momenta, A£J oc k^ . 

Figures 01a) and^b) respectively show the calculated 
electron spin coherence time and measured mobilitjsS as 
a function of temperature for bulk InAs. At 300 K, T2 is 
20.6 ps for /i = 2.89 x 10"* cm^/Vs, and it is one order in 
magnitude shorter than the electron spin coherence time 
of GaAs. This is consistent with the recent measurements 
of n-type InAs at room temperature by Boggess et ali^. 
The measured mobility of the sample in Ref. [13 is 2.6 x 
10^ cm^/Vs and the electron density is approximately 
2 X 10^^ cm~^. The measured Ti is 19 ± 4 ps which is in 
excellent agreement with the calculated Ti of 21 ps. 

Shown in Figs. 0Jc) and Eld) are the calculated elec- 
tron spin coherence time and measured mobility^ as a 
function of temperature for bulk InSb, respectively. The 
calculated T2 is 11.6 ps at 300 K, and the measured mo- 
bility is 1.37 X lO'' cm^/Vs. The electron spin coherence 
time of InSb is shorter than that of InAs at 300 K due to 
the larger conduction spin splitting in InSb [see Figs.|2ta) 
and|3{b)]. Note that Ti and T2 scale approximately as 
/i"^. The band structures of InAs and InSb in our model 
are shown in Figs. |3a)-(b). The ratio of the spin-orbit 
coupling to the band gap i?inAs = 1-07 for InAs and i?inSb 
= 5.04 for InSb. 

The calculated electron spin coherence time as a func- 
tion of temperature for bulk InP is shown in Fig. IHI^a) 
and the temperature-dependence mobility^*' from exper- 
iments is depicted in Fig. EJb). The calculated electron 
spin coherence time T2 is 0.49 ns at 300 K for /i = 
3.90 X 10^ cmVVs with i?i„p = 0.08 [see Fig. 0. The 
longer electron spin coherence times in InP come from 
the smaller conduction spin splittings [see Figs.|3|(c) and 
(d)]. 



2. 11- VI SEMICOND UCTORS 

Recently II- VI zincblende semiconductors have also 
gained attention due to their long electron spin coher- 
ence times and transport properties. It has been exper- 
imentally demonstrated by Malajovich et al. that spin 
information can be protected by transport from a region 
of high spin decoherence (GaAs) to a region (ZnSe) of 
low spin decoherence^^. Here we calculate the electron 
spin coherence times of bulk II- VI semiconductors and 
we will look at the electron spin coherence times of II- VI 
quantum-well semiconductors in the next section. 

Figures Eta) and |HIb) show respectively the electron 
spin coherence time T2 and measured mobility as a func- 
tion of temperature for bulk CdSe. Due to the relatively 
small spin-orbit coupling (see Table IVH) . the calculated 
electron spin coherence time is very long and T2 — 96 
ns at 300 K with fi = 600 cm^/Vs. For bulk ZnSe, the 
calculated electron spin coherence time T2 = 25 ns at 300 
K with /i = 388 cm^/Vs. The calculated electron spin 
coherence time T2 and measured mobility as a function of 
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temperature for ZnSe are respectively shown in Figs.|HIc) 
andlH^d). The electronic band structures for CdSe and 
ZnSe are displayed in Figs.|nta)-(b), and the spin-orbit 
coupling to band-gap ratio for CdSe and ZnSe is respec- 
tively i?cdSc — 0.15 and i?znSc — 0.26. We approximate 
both systems in our model as zincblende lattices. 

For bulk semiconductors the relevant electronic states 
for electron spin decoherence are the states of the two 
lowest conduction bands near the bulk band edge, and 
thus perturbative expressions22i^ for ri^;(l, E) for these 
bulk semiconductors [fl^^^{l,E) oc E^] are nearly iden- 
tical to those obtained from a full fourteen-band calcu- 
lation (within numerical accuracy). The energy depen- 
dence of hl33{l,E) for the III-V and II- VI bulk semi- 
conductors considered above is summarized in Fig. 1101 
and the cubic dependence on energy of ^1^^{l,E) [i.e., 
the slopes of the curves are approximately 3 and thus 
^£c33(l, k) ex k'^] for bulk semiconductors near zone cen- 
ter is confirmed. The calculated spin coherence times for 
bulk semiconductors using the perturbative expression 
obtained by D'yakonov and Perel™ coincide with those 
calculated using the fourteen-band non-perturbative the- 
ory. The agreement between calculated and measured 
spin coherence times (T2's) indicates that the conduction- 
band electron spin splitting of bulk zincblende semi- 
conductors is well described by the fourteen-band non- 
perturbative theory. 



SPIN COHERENCE TIMES IN QUANTUM 
WELLS AND SUPERLATTICES 



Electron spin dynamics in III-V and II- VI quantum 
wells have been extensively investigated experimentally 
because of their well-known electronic and optical proper- 
tics. Ultrafast all-optical spin polarization switches using 
multiple quantum wells have been demonstrated exper- 
imentally by Nishikawa et alm^, and the result suggests 
that these all-optical spin-dependent switches can be 
used for future ultrafast optical communication systems 
that require very high switching speeds. The switching 
performance of these all-optical spin-dependent switches 
hinges on the electron spin decoherence, and therefore, 
the calculations of electron spin coherence times are im- 
portant. 



1 . Ill- V LA YERED SEMICOND UCTORS 

Terauchi et ali^ have performed measurements on 
GaAs/AlGaAs multiple quantum wells to study the mo- 
bility dependent electron spin coherence times. The 
band-edge diagram and electronic band structure of a 
75-A GaAs/100-A Alo.4Gao.6As quantum well are shown 
in Figs. IllT a^-fb). and the calculated band-gap energy 
is 1489 meV (0.83 /im). The lowest conduction-band 
spin splitting as a function of crystal momentum k||[110] 



and excess energy for this quantum well is plotted in 
Figs. IT^ a'l-l'b') (along with other structures considered 
later). Figure ITST a) shows the calculated electron spin 
coherence times Ti 's as a function of mobility assuming 
optical phonon (solid line) and neutral impurity (dashed 
line) scattering as the dominant process determining the 
mobility in the quantum well at 300 K. In the figure, the 
results of the experimental measurements are represented 
by closed circles, and, for comparison, the theoretical re- 
sults of the DK theory based on perturbative approach 
are also shown. Our results agree with experiment if one 
assumes a shift from optical phonon to neutral impurity 
scattering as the mobility drops — this is the origin of the 
unusual experimental dependence of Ti on the mobility. 
The calculated electron spin coherence times with opti- 
cal phonon scattering are roughly two times larger than 
those with neutral impurity scattering, and this is di- 
rectly related to the energy-dependence of the transport 
times for different momentum scattering mechanisms (see 
Table ini). 

Note that the experimental results have been adjusted 
from Ref . llSl for the authors defined an effective spin flip 
time for a single spin, t^ = 2Ti , and plotted their results 
for Tg . The discrepcncy in the magnitude of Ti 's between 
the DK calculations and measurements is about a factor 
of 4 (see Ref. l20j) . The electron spin coherence times Ti's 
as a function of average transport time Tp and electron 
density n are displayed in Figs. ll3r b') and ll3r c'). respec- 
tively. Open squares and triangles respectively represent 
the calculated Ti's for optical phonon and neutral im- 
purity scattering mechanisms. The relationship between 
the measured electron mobility and density is depicted 
inFig.inid). 

Shown in Fig. E| are the calculated electron spin co- 
herence times Ti 's and the corresponding spin relaxation 
rates T^ 's as a function of electron confinement energy 
Eie for (30 - 100)-A GaAs/Alo.4Gao.6As quantum wells 
at 300 K assuming optical phonon (solid line) and neutral 
impurity (dashed line) scattering. The results of experi- 
mental measurements-^ are represented by closed circles 
and the calculations of DK theory are represented by 
dot-dashed line. Our numerical results are in excellent 
agreement with the experimental measurements both in 
the absolute magnitude and trend. On the other hand, 
the calculations of the DK theory differ approximately by 
a factor of 4 and the predicted trend is stronger than the 
observed ona^. The electron spin coherence times and 
electron spin relaxation rates as a function of quantum- 
well width are illustrated in Fig. [HI 

Another set of comprehensive measurements of the 
electron spin coherence time Ti for (25 — 300)-A 
GaAs/150-A Alo.35Gao.65As quantum wells was recently 
carried out by Malinowski et al}'^, and the results of the 
experiments are shown in Fig. E| The calculated Ti's 
as a function of electron confinement energy for electron 
density n = 1 x 10^^ and an estimated constant electron 
mobility /x = 4600 cm^/Vs are in reasonable agreement 
with the experimental measurements, in particular for 
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quantum wells with relatively large electron confinement 
energy Ei^. The slower measured electron spin relax- 
ation rates T^ 's are most likely associated with the non- 
uniform electron mobility in different thickness quantum 
wells. As the quantum- well width increases, the calcu- 
lated electron spin coherence times approach the values 
of bulk GaAs. The calculated spin coherence times of 
bulk GaAs at room temperature with n = 1 x 10^^ and /i 
= 4600 cm^/Vs for optical phonon and neutral impurity 
scatterings are approximately 29 ps and 27 ps, respec- 
tively. Figure 1161 shows the temperature dependence of 
Ti for three different width quantum wells (60-A, 100- 
A, and 200-A) with the same electron density and esti- 
mated mobility as in Figs. 1151 Given the uncertainty in 
the electron mobility, the calculated Ti's compared with 
the results of experiments are quite reasonable. 

The electron spin coherence times Ti's and the cor- 
responding spin relaxation rates Tf^'s of (25 — 100)-A 
Ino.53Gao.47As/97-A InP quantum wellsii as a function 
confinement energy Ei^ and quantum-well width LwoU 
at 300 K for electron densities n = 6 x 10^^ cm~'^ and 
n = 5 X 10^^ cm~^ are shown in Fig. El Since the ex- 
perimental electron mobility is not available, a typical 
value of 6700 cm^/Vs is used in our calculations. The 
confinement-energy dependence of Ti's from our theo- 
retical calculations agrees with those of the experimen- 
tal measurements. On other hand, as in the case of 
GaAs/AlGaAs quantum wells, the DK theory gives a 
stronger confinement-energy dependence of the electron 
spin coherence time than that of the experimental results. 
Better agreement between our theorectical calculations 
and experimental results is obtained if the electron den- 
sity is degenerate. 

The band-edge diagram and electronic band structure 
for a 70-A Ino.53Gao.47As/97-A InP quantum well are 
illustrated in Figs. I18f a)-(b) with band-gap energy of 
976 meV (1.27 fim). The conduction-band spin splitting 
of the lowest states as a function of crystal momentum 
k||[110] and excess energy is shown in Figs. Efa)-(b). 
The measured Ti for this quantum well by Tackeuchi et 
alm^ is 2.6 ps, and calculated Ti is 4.0 ps (6.5 ps) using 
an estimated electron mobility /i — 6700 cm^/Vs and a 
value of degenerate electron density n — 3.0 x 10^^ cm~^ 
assuming neutral impurity (optical phonon) scatteringSS,. 
InGaAs/InP quantum wells are a so called non-common- 
atom system^, which gives rise to interface bonding 
asymmetry or native interface asymmetry (NIA) . Native 
interface asymmetry causes additional electron spin split- 
tings, and consequently it leads to faster electron spin 
decoherence in the system. We do not include NIA con- 
tributions in our calculations for these wells. 

We also present the electron spin coherence times 
for InGaAsP quaternary quantum wells with different 
compositions. The calculated Ti for neutral impurity 
scattering at 300 K with electron density n = 1.3 x 10^^ 
cm~'^ and electron mobility fj, = 800 cm'^/Vs for a 55-A 
Ino.53Gao.47As/lOO-A Ino.87Gao.13Aso.29Po. 71, 65-A 
Ino.53Gao.47As/lOO-A Ino.87Gao.13Aso.29Po. 71, 90-A 



In0.53Ga0.47As0.93P0.07/75-A Ino.87Gao.13Aso.29Po.71, 
and 95-A In0.53Ga0.47As0. 93P0.0 7/75- A 

Ino.87Gao.13Aso.29Po. 71 is listed in Table IVIIII The 
calculated electron spin coherence times of these quan- 
tum wells are comparible to those of InGaAs/InP 
quantum wells [c.f., Fig. I17| . These results are consis- 
tent with the experimental measurements by Hyland 
et al. in Ref. 13 (see Table NTU^ . The band-edge 
diagram and electronic band structure for a 95-A 
In0.53Ga0.47As0.93P0.07/75-A Ino.87Gao.13Aso.29Po.71 
quantum well, and the band-gap energy of 830 meV 
(1.49 ^m) are displayed in Figs. Ela)-(b). The lowest 
conduction-band spin splitting as a function of crystal 
momentum k||[110] and excess energy is depicted in 
Figs.[Ilc)-(d). 

The confinement-energy dependence of the elec- 
tron spin coherence times at room temperature for 
GaSb/AlSb and GaAsSb/AlSb quantum wells with elec- 
tron density n = 2 x 10^^ cm~"^ and estimated mobility 
/i = 2000 cm^/Vs are respectively shown in Figs. I20r a') 
and l20f b). The Ti's for both structures have very similar 
confinement-energy dependence, and the coherence times 
for the GaAsSb/AlSb quantum wells are slightly longer 
than those of the GaSb/AlSb quantum wells. The longer 
electron spin coherence times in the GaSb/AlSb quan- 
tum wells are attributed to the smaller conduction-band 
spin splittings due to the 19% indium substitution in the 
well-acting layers. 

The band-edge diagram and electronic band struc- 
ture for a 80- A GaSb /80-A AlSb QW and a 51- A 
Gao.19Aso.8iSb /80-A AlSb QW are shown in Figs.lHTa)- 
(b) and|21Ia)-(b), respectively. For the former, the band- 
gap energy is 829 meV (1.50 /J.m), while for the latter 
the band-gap energy is 795 meV (1.56 /.tm). The spin 
splitting of the lowest conduction band as a function of 
crystal momentum k| | [110] and excess energy is plotted in 
Figs.Efa)-(b). Hall et alm^ have experimentally investi- 
gated the electron spin dynamics for both quantum wells, 
and measured Ti at 295 K is 0.52 ps for the GaSb/AlSb 
QW with optical excited electron density n — 2.8 x 10^® 
cm~^, and for the GaAsSb/AlSb QW, the measured Ti 
at 295 K is 0.42 ps with optical excited electron density 
n = 3.4 X 10"'^* cm""^. The calculated electron spin coher- 
ence time Ti for the GaSb/AlSb QW is 0.53 ps (1.4 ps) 
assuming neutral impurity (optical phonon) scattering; 
and the calculated Ti for the GaAsSb/AlSb QW is 0.87 
ps (1.4 ps) assuming neutral impurity (optical phonon) 
scatteringSS. The estimated electron mobility is /i = 2000 
cm^/Vs for both quantum wells. In this case, the agree- 
ment between theory and experiment is quite good. It 
should be noted that the optical excited electron densities 
in these experiments are highly degenerate. 

Due to the relatively large electron spin precession vec- 
tor of GaSb, the electron spin coherence times in GaSb 
quantum wells are expected to be relative short, for ex- 
ample, compared to the electron spin coherence times of 
GaAs quantum wells [c.f.. Figs. UlTa) and EUTa)] . Fig- 
ure I23r a') displays the dependence of the electron spin 
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coherence times on temperature for a 80-A GaSb/300- 
A Alo.4Gao.6Sb QW. The experimental measurements 
of electron mobilities as a function of temperature are 
shown in Fig. l23r b). At room temperature, the expected 
electron spin coherence time Ti = 1.2 ps for /i — 1300 
cm^/Vs. The electron spin coherence times are propor- 
tional to the electron mobility but inversely proportional 
to the temperature. This example illustrates that in or- 
der to determine the dominant electron spin relaxation 
mechanism in a system, quantitative calculations of the 
electron spin coherence times are necessary for the elec- 
tron spin coherence times depend on the electron mobil- 
ity, density, and temperature. 

Next we calculate the electron spin coherence times 
at room temperature for a thin-layer type-II superlat- 
tice, which is a non-common-atom structure. The cal- 
culated Ti's for a 21.2A InAs/36.6A GaSb superlattice 
are respectively 0.77 ps and 1.6 ps for optical phonon 
and neutral scattering with electron density n — 2x 10^^ 
cm"'^ and mobility fi = 3000 cni'^/Vs. The band-edge 
diagram and electronic band structure are displayed in 
Figs. |23Ia)~(b) with band-gap energy of 303 meV (4.09 
/im). The lowest conduction-band spin splitting as a 
function of crystal momentum k||[110] and excess en- 
ergy for the superlattice is displayed in Figs. EJa)-(b). 
The effects due to native interface asymmetry are not in- 
cluded in these calculations, and it has been shown by 
Olesberg et alm^ that native interface asymmetry plays 
an important role in electron spin decoherence in this 
particular structure (for details see Ref. l63). 

Another non-common-atom type-II system of inter- 
est in spintronic devices is InAs/AlSb quantum wells in 
which the conduction-band offset AEc is approximately 
one order in magnitude larger than the valence-band 
offset AEv [Fig. I25r a')]. The band-edge diagram and 
electronic band structure for a 75-A InAs/200-A AlSb 
quantum well are illustrated in Figs. l^ aVfb). and the 
band-gap energy is approximately 300 meV (4.13 /im). 
The confinement-energy dependence of the electron spin 
coherence times at room temperature for (25 — 100)- 
A InAs/200-A AlSb quantum wells with electron den- 
sity n — 2 X 10^^ cm^"^ and an estimated mobility fi = 
2000 cm^/Vs is shown in Fig. l26f a). and the correspond- 
ing electron spin relaxation rate Tj~ as a function of 
quantum- well width iwcU is displayed in Fig. EHb) . The 
range of tunability of the electron spin coherence time in 
this system is quite large, over two orders of magnitude. 
For a wide quantum well, e.g., -Lwcii = 300A Ti is ap- 
proximately 70 ps; in contrast, Ti is approximately 0.2 
ps for a narrow quantum well of width Lwoii — 30 A. 

The electron spin coherence times as a function of elec- 
tron mobility and density for a 75-A InAs/200-A AlSb 
quantum well are shown in Figs. l26f c)-fd'). As expected, 
for non-degenerate electron density, Ti oc /i^^. It should 
be noted that, for systems with non-degenerate carrier 
density, the carriers in the conduction band occupy rel- 
atively low energy states (-Boccupiod < k^T). Within 
this energy range, the electron spin coherence time de- 



creases with increasing carrier density because the effec- 
tive magnetic field increases with increasing energy as 
the higher energy states are occupied [see Fig. EHd)]. 
For degenerate systems, the electron spin coherence time 
may decrease or increase with increasing carrier density 
depending on where the Fermi level lies because the ef- 
fective magnetic field oscillates as the energy increases 
[see Figs. |27|;a)-(d)]. 



2. SPIN-LED STRUCTURES 

Another technique used to detect the spin polariza- 
tion of carriers in a semiconductor spintronic device is 
to measure the degree of spin polarization of the electro- 
luminescence (EL) in a light-emiting diode (LED) fabri- 
cated within the deviceMiMi^. This technique provides 
a quantitative and model-independent measurement of 
electrical spin injection^. A crucial factor in determin- 
ing the spin polarization is the carrier spin relaxation 
after spin injection; therefore, it is important to know 
the spin coherence times (both Ti and T2) of the carri- 
ers in order to optimize the design of the device. In the 
following, we will calculate the spin coherence time for 
some common spin-LED devices used to investigate the 
efficiency of electrical spin injection. 

The electronic band structure and band-edge diagram 
for a 100-A Ino.13Gao.87As/2OO-A GaAs QW^^ are shown 
in Figs. lTTT aVfb). and the band-gap energy is 1326 meV 
(0.94 /im). The lowest electron spin splitting of the low- 
est states as a function of crystal momentum k||[110] 
and excess energy for the quantum well is displayed in 
Figs. Elc)-(d). Figure OST al shows the mobility de- 
pendence of the spin coherence times at room temper- 
ature. The Ti's for optical phonon scattering differ from 
those for neutral impurity scattering by a factor of 2, 
and the inverse relationship between the electron spin 
coherence time and mobility is the signature of DP spin 
relaxation mechanism for systems with a non-degenerate 
carrier density [see Eq. (|82a|l ]. The electron spin coher- 
ence times as a function of carrier density are illustrated 
in Fig. l28r b). and the non-linear behavior of the electron 
spin coherence times are due to T^^ c>c ^1i{C',E) [see 
Figs.|S7i;a)-(d)]. 

The band-edge diagram and electronic band struc- 
ture for a 40-A Ino.20Gao.80As/lOO-A GaAs/40-A 
Ino.2oGao.8oAs/200-A GaAs double quantum well^^ at 
room temperature are depicted in Figs. EW a'l-fb) with 
band-gap energy of 1322 meV (0.94 /im). The lowest 
conduction-band spin splitting as a function of k||[110] 
and excess energy for the quantum well is plotted in 
Figs. ll2f c')-('d). The mobility and carrier-density depen- 
dence of the spin coherence times at room temperature 
are respectively shown in Figs. EHIc) and EST d). Simi- 
lar dependences of the electron spin coherence times on 
mobility and carrier density are found. 

Next we would like to compare the electron spin co- 
herence times of two spin-LED structures with differ- 
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ent aluminum concentration in the barriers for otherwise 
identical quantum wells. The band-edge diagrams and 
electronic band structures for these two structures are 
shown in Figs. I30f a)-fb) and|^^a)-(b). The band-gap 
energy for the spin-LED structure shown in Figs. l30f a)- 
(b) is 1432 meV (0.87 /im), and the band-gap energy for 
the spin-LED structure displayed in Figs. I3ir a)-fb) is 
1445 meV (0.86 /im). The conduction-band spin splitting 
of the lowest states as a function of crystal momentum 
k||[110] and excess energy for these two spin-LED struc- 
tures is illustrated in Figs.EJc)-(d). The mobility and 
carrier-density dependence of the electron spin coherence 
times for a 150-A GaAs/100-A Alo.03Gao.97 As QW^^^ at 
room temperature is plotted respectively in Figs. \^ a,) 
and l32r b^. while the mobility and carrier-density depen- 
dence of the electron spin coherence times for a 150-A 
GaAs/100-A Alo.4Gao.6As QW at room temperature is 
respectively shown in Figs.|22Ic) and l32f d). The electron 
spin coherence times for the quantum well with 3% alu- 
minum concentration in the barriers are approximately 
a factor of 5 longer than those of the quantum well with 
40% aluminum concentration in the barriers. As the alu- 
minum concentration increases, the penetration of the 
electron wavefunctions into the barriers decreases due 
to the increase of potential energy in the barriers. As 
a result, the conduction spin splittings increase as the 
confinement energy increases [see Figs.|^Jc) and l37r d')]. 
Consequently the electron spin coherence times decrease 
as the aluminum concentration increases. 

Next we consider a thin 50-A GaAs quantum well with 
100-A Gao.51Ino.49P as barriers. Shown in Figs.|23a)-(c) 
are the dependence of the electron spin coherence times 
on the electron mobility, average electron momentum 
scattering time, and temperature. Although the tem- 
perature dependence of the electron mobility is similar 
to that of the GaSb quantum well considered above [see 
Figs. HSId) and lSST d)]. the dependence of the electron 
spin coherence times on temperature is quite different 
[c.f., Figs.[53Jc) and l33r c)]. In this case the electron spin 
coherence times are proportional to the electron mobil- 
ity at high temperature and inversely proportional to the 
electron mobility at low temperature. 



3. II- VI LA YERED SEMICOND UCTORS 

Recently it has been demonstrated experimentally by 
Kikkawa et al^ that II- VI semiconductor quantum wells 
have remarkably long electron spin coherence times, of 
the order of nanoseconds, which persist to room temper- 
ature. Here we look at the electron spin decoherence of 
two II- VI semiconductor quantum wells: CdSe/ZnSe and 
CdZnSc/ZnSc quantum wells. The electron spin coher- 
ence times as a function of mobility, average momentum 
scattering time, and temperature for a 105-A CdSe/500- 
A ZnSe quantum well-^ are plotted in Figs. \'64( a,)-(c). 
The electron mobilities are extrapolated from a mea- 
sured value of 6800 cm^/Vs at 4.2 K assuming neutral 



impurity or optical phonon scattering, and the results 
are displayed in Fig. ISlT d'). The calculated electron spin 
coherence times are proportional to the temperature, but 
they are inversely proportional to the electron mobility 
because the extrapolated electron mobilities are inversely 
proportional to the temperature. The calculated electron 
spin coherence times at room temperature are approxi- 
mately 4.7 ns and 2.7 ns for optical phonon and neutral 
impurity scattering, respectively. The calculated elec- 
tron spin coherence times are comparable to those of the 
experimental measurements reported in Ref. |j. 

Shown in Figs. ISCT al-fbl are the band-edge di- 
agram and electronic band structure for a 105-A 
Cdo.24Zno.76Se/5OO-A ZnSe quantum wellSi with band- 
gap energy of 2425 meV (0.51 /im). Figures l^ al-fc) 
illustrate the dependence of the electron spin coherence 
times on electron mobility, average momentum scattering 
time, and temperature. The temperature dependence of 
the electron mobilities are shown in Fig. l36r d). which is 
obtained by extrapolating the electron mobility from a 
measured value of 8800 cm^/Vs at 4.2 K assuming neu- 
tral impurity and optical phonon scattering. The depen- 
dence of the electron spin coherence times on mobility, 
average momentum scattering time, and temperature is 
very similar to that of the CdSe quantum well considered 
above. 



V. VALIDITY OF DK APPROXIMATIONS 

The origin of the discrepancies between the DK theo- 
retical and the experimental results can be traced back 
to the approximations made in the derivation of Ti in 
the DK theory. To summarize, the expression for the 
electron spin coherence time Ti in the DK theory for 
(OOl)-quantum wells is obtained as follows (see Sec llIIfl . 
First, it is assumed that the penetration of the electronic 
wavefunctions into the barriers in the quantum well is 
negligible, hence Eic ^ Ai?c, where Eic is the electron 
confinement energy of the first quantum-well state and 
AEc is the conduction band offset. Second, the per- 
turbative expressions^ Clli{l,E) ex E{4:Eie - E)"^ and 
0^3(1, i?) oc E^ are used. Third, the energies of rele- 
vant states are assumed to be ii^ ^ Eic, and thus (i) the 
contribution from ^^si^iE) is ignored, and (ii) it is as- 
sumed that ^^1(1, E) (X E. It is not generally recognized 
that the conditions /csT <C i?ic ^ ASc are quite restric- 
tive and are difficult to satisfy at room temperature. It 
should be noted that in this approximation only the first 
quantum- well state is occupied, and hence it can be only 
applied to quantum wells with large energy separation 
between the first and the second quantized states. 

Under the DK assumptions, the electron spin coher- 
ence time Ti is inversely proportional to the mobility 
independent of the dominant scattering mechanism. In 
addition, Ti is inversely proportional to the confinement 
energy Ef^. These trends are not supported by recent 
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experimental measurementg"'^^ illustrated in Figs. llSf a) 
and I14r a') . In both cases the experimental trends are 
weaker than the predicted theoretical ones. The weaker 
dependence of Ti on Eic in our theory versus DK theory 
is due to wavefunction penetration into the barriers and 
non-perturbative effects (the spin splitting for damped 
waves in the barrier is of the opposite sign as that for 
propagating waves in the well) . 

The validity of DK approximations can be further ana- 
lyzed quantitatively by comparing the energy dependence 
of fl'^i{l,E) and fl^^{l,E) for these material systems. 
In DK theory, i^^i(l,^} is a linear function of energy E 
(see SecHniand Ref.'si'); however, Figs.Eta) andEZ^c) 
show that for quantum wells Q,^i{l, E) is only linear for 
a small energy range {E < k-gT) above the band edge be- 
fore it begins to deviate [c.f., the short dashed line based 
on the DK theory and the solid line in Fig.|23a)]. More 
energetic states than this certainly contribute to the elec- 
tron spin decoherence at room temperature. The wider 
the well the lower the energy where Vi^^{l,E) deviates 
from linear behavior [c.f., the solid line in Figs. lSTf a') and 
ISTT c)]. as it approaches a bulk-hke E'^ behavior [Fig.[TU|. 
r2^3(l,i?) for the material systems considered above is 
shown in Figs. EZI^b) and I37f d). and for a large width 
well is comparable in magnitude to ^^i(l, E). 

The applicability of the perturbative approach for 
calculations of electron spin coherence times in two- 
dimensional systems can be explored by carefully examin- 
ing the energy dependence of Q,xi{\, E) and ^^^{l^ E) for 
various structures. As the wells become narrower, even 
the perturbative expressions for Vlxi{^,E) and D,x3{l, E) 
break down. Figures E^a) and OTTb) show VlI^{1,E) 
and 51^3(1, E) (long dashed lines) for a thin-layer type-II 
InAs/GaSb superlattice, indicating very different behav- 
ior from the other structures, poorly reproduced by even 
the general forms of the perturbative expression. 



important role of momentum scattering mechanisms in 
determining the electron spin coherence times by study- 
ing the mobility dependence of Tx for various momentum 
scattering mechanisms. 
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VI. CONCLUDING REMARKS 



We have presented a quantitatively accurate non- 
perturbative nanostructure theory for electron spin de- 
coherence in bulk and quantum- well zincblende semicon- 
ductors based on a fourteen-band model. The calcu- 
lated electron spin coherence times in bulk and layered 
zincblende semiconductors are in agreement with exper- 
imental measurements, indicating the importance of ac- 
curate band-structure calculations for zincblende type 
nanostructures. We have explicitly identified the im- 
portant assumptions of the DK theory regarding both 
quantum- well electronic structure and carrier dynamics, 
and the origins of the discrepancies between the results 
of the DK theory and a full fourteen-band theory. The 
validity and limitation of the perturbative approach are 
investigated numerically through specific examples, and 
the results show that the perturbative approach breaks 
down for systems with narrow wells. We have shown the 
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TABLE I: The ratio rj(^^^ of the transport time ri(n, fc) and the momentum relaxation time 
Ti{n,k), and the power law of the energy-dependent transport time [ri(n, fc) oc E^^Ck)] for 
ionized impurity, optical phonon, and neutral impurity scattering in bulk semiconductors'^". 
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TABLE IL The ratio r;^^-^) of the transport time ti{C,K\\^,Kz) and the momentum relaxation time 
ti{CtK\\,Kz) for different angular index Z, and the power law of the energy-dependent transport 
time \ti{C,K\\,Kz) oc ^^^(K)] for various scattering mechanisms in quantum wells or superlatticea^. 
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TABLE IIL Bulk basis states |n,(T, 0) used in the fourteen-band model for the Fe 
bulk band edge and the bonding and antibonding r7 and Fg bulk band edges. 
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TABLE V: Parameters for fourteen-band k-p band structures for Alj;Gai_a;Sb, A^Gai-^^As, Ina;Gai_a;As, GaAs2:Sbi_a;, and Gaa;Ini„a;P at room tempera- 
ture. Eg is the fundamental gap (between s and p (Fg) states) , Ec is the zone-center energy difference between p (Fg) and p (F7) , and A and Ac are the spin- 
orbit spUtting of the p and p states, respectively, Co is the lattice constant, cn and C12 are the elastic constants, and Va and Vb arc the deformation potentials. 



ALGai-:,Sb 



ALGai-:cAs 



Inj:Gai-iAs 



GaAsa;Sbi 



Gaxlni-xP 



0.724 + 1.105x-h0.471x^ 

0.673a; + 0.758(1 - x) 

3.50a; + 3.19(1 -a;) 

0.0a; + 0.21(1 - a;) 
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12.0a; + 15.1(1 -a;) 

6.136a; + 6.096(1 ~ x) 

8.769a; + 8.834(1 - x) 

4.341a; +4.023(1 - x) 

-5.9a; -7.2(1 -a;) 

-1.4x- 1.8(1 -x) 



0.724 - 0.506a; + 1.204a;^ 1.351 + 0.647a; + 0.7859a;' 
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A (eV) 
Ec (eV) 
Ac (eV) 
EPo (eV) 
Ep, (eV) 
Eq (eV) 

fflo (A) 
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C12 (10" dyn/cm^) 

Va (eV) 

14 (eV) 



1.422 + 1.198X + 0.433a;' 

0.30X + 0.341(1 -x) 

4.48a,- + 4.61(1 - x) 

0.0a; + 0.16(1 -x) 

23.5x + 27.0(1 - x) 

0.03X + 5.9(1 -x) 
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12.02x + 12.3(1 - x) 

5.70X + 5.71(1 -x) 

-2.6x- 9.8(1 -x) 

-1.5x- 2.0(1 -x) 



1.422 - 1.484x + 0.418x' 
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21.6X + 27.0(1 -x) 
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15.8X + 12.8(1 -x) 
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-1.8x- 2.0(1 -x) 
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4.77X + 3.19(1 -x) 

0.16X + 0.21(1 -x) 

27.0X + 27.5(1 - x) 

5.90X + 11.5(1 -x) 

12.8X + 15.1(1 -x) 

5.653X + 6.096(1 -x 

12.3x + 8.834(1 -x) 

5.71X + 4.023(1 -x) 

-9.8x- 7.2(1 -x) 

-2.0x- 1.8(1 -x) 



0.080X + 0.108(1 -x) 
5.24X + 4.80(1 -x) 
0.09X + 0.07(1 -x) 
22.0X + 20.2(1 -x) 

3.0X + 3.7(1 -x) 

22.8X + 15.9(1 -x) 

5.447X + 5.860(1 - x) 

14.388X + 10.11(1 -x) 

7.144X + 5.61(1 -x) 

-9.3x- 6.4(1 -x) 

-1.8x- 2.0(1 -x) 
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TABLE VI: Parameters for fourteen-band k • p band structures for GaP^jAsi-a;, InPa;Asi_a;, In2:Gai_a;Sb, and Cda;Zni„a;Se at room temperature. Eg is 
the fundamental gap (between s and p(Xs) states), Ec is the zone-center energy difference between ^(rs) and p{T-/), and A and Ac are the spin-orbit 
sphtting of the p and p states, respectively, Uo is the lattice constant, cn and C12 are the elastic constants, and Va and Vb are the deformation potentials. 

GaPaiAsi-j, InPjjAsi-a: luxGai-xSh CdaiZni-j^Se 

Eg (eV) 1.422 + 1.172a; -hO.igOa;^ 0.356 + 0.675a; + 0.320a;'' 0.723 - 0.926a: + 0.371a;^ 1.670 - 0.673a; + 0.302a;^ 

A (cV) O.OSOx + 0.341(1 - a;) 0.108a; -h 0.380(1 - a;) 0.850a; -h 0.758(1 - a;) 0.390a; + 0.417(1 - a;) 

Ec (cY) 5.24x -H 4.61(1 - a;) 4.63a; -h 4.39(1 - a;) 3.14a; -F 3.19(1 - a;) 6.50a; -h 7.17(1 - a;) 

Ac (eV) 0.09a; + 0.16(1 -a,-) 0.07a; + 0.24(1 - a;) 0.39a; + 0.21(1 - a;) 0.24a; + 0.08(1 - a;) 

Spo (eV) 22.0a; + 27.0(1 - a,-) 20.2a; + 21.6(1 - a;) 24.0a; + 27.5(1 - a;) 23.1a; + 19.9(1 - a;) 

Ep, (eV) 3.0a; + 5.9(1 - a;) 3.7a; + 5.2(1 - a;) 7.5a; + 11.5(1 - a;) 0.05a; + 0.2(1 - a;) 

Eq (eV) 22.8a; + 12.8(1 - x) 15.9a; + 15.8(1 - a;) 15.5a; + 15.1(1 - a;) 13.0a; + 12.9(1 - a;) 

tto (A) 5.447a; + 5.653(1 - a;) 5.860a; + 6.058(1 - x) 6.469a; + 6.096(1 - x) 6.077a; + 5.668(1 - a;) 

cii (10" dyn/cm^) 14.388a; + 12.3(1 - a;) lO.llx + 8.329(1 - a;) 6.918a; + 8.834(1 - a;) 7.49x + 9.407(1 - a;) 

C12 (10" dyn/cm^) 7.144a; + 5.71(1 - a;) 5.61a; + 4.526(1 - a;) 3.788a; + 4.023(1 - a;) 4.609a; + 5.625(1 - a;) 

K (eV) -9.3a; - 9.8(1 - a;) -6.4a; - 5.8(1 - a;) -7.7a; - 7.2(1 - x) -2.1a; - 6.2(1 - a;) 

Vt (eV) -1.8a; -2.0(1 -a;) -2.0a; - 1.8(1 - a;) -2.0x - 1.8(1 - x) -l.lx - 1.4(1 - a;) 

Epg = 2mc^PQ, Ep^ = 2me^Pi, and Eq = 2m^^Q2 



to 



TABLE VII: The conduction-band offset at 300 K for various structures 



Structure 



75-A GaAs/100-A Alo.4Gao.6As QW 

150-A GaAs/100-A Alo.03Gao.g7 As QW 

150-A GaAs/100-A Alo.4Gao.6As QW 

21.2-A InAs/36.6-A GaSb SL 

75-A InAs/200-A AlSb QW 

70-A Ino.53Gao.47As/97-A InP QW 

100-A Ino.13Gao.87As/2OO-A GaAs QW 

95-A In0.57Ga0.43As0.93P0.07/75-A Ino.87Gao.13Aso.29Po.71 QW 

40-A Ino.20Gao.80As/lOO-A GaAs/40-A Ino.20Gao.80As/2OO-A GaAs QW 

80-A GaSb /80-A AlSb QW 

51-A Gao.019Aso.8iSb /80-A AlSb QW 

105-A Cdo.24Zno.76Se/75-A ZnSe QW 



AEc (meV) 



AE, (moV) 



370 

22 
360 
981 
2188 
261 

78 
161 
121 
1212 
1511 
114 



4 



4 
49 
106 


62 

2 

99 
35 
105 
82 



AEc is the conduction- band offset, and AEs is the spin splitting energy due to lattice-mismatched-induced strain. 
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FIG. 1: The fourteen-band electronic band struc- 
ture of bulk III-V semiconductors near the T point 
of the BriUiouin zone at 300K. (a) GaAs. (b) GaSb. 
The fourteen bands considered are six p states, two 
s states, and six p states (s indicates antibonding 
s states and p indicates antibonding p states). The 
s — p, s — p, and p — p interband interactions are 
characterized by the corresponding interband mo- 
mentum matrix elements Pq: Piy ^^id Q, respectively. 
The energy parameters Eg, Ec, A, and Ac are the 
energy difference of s (Fg) - p (Fs), p (Fy) - p (Fs), 
Pi^s) - PiTr), and ^(Fs) - ^(Ft), respectively. 



10^ 


■(a) 


"1 


^in^ 


r \. 


-1 


"^ 5 


^^ 


: 


S 10 

h."in' 


1 1 

/ 


- 


in^ 


r ^X X 


■ T 


10^ 


N 




^ "^^ 


_^ 


.„i 


-, , 1 1 


, - 



r(K) 

10 100 300 



in^ 


^c) 


' 1 1 


Iio' 


^ 


^ - 


io' 


r ^ 


^^^^1 


.„o 




, 1 1 



10 




1 00 200 

r(K) 



100 200 300 

r(K) 



FIG. 2: Temperature dependence of mobility and 
coherence time in bulk III-V semiconductors, (a) 
T2 as a function of temperature for GaAs with n = 
1.0 X IQis cm~3 (soii(j line^O) and n = 1.0 x 10" 
cm^'^ (dashed line), and squares represent the re- 
sults of experiments^, (b) fi vs T for GaAs with 
n = 1.0 X 10^^ cm-3 (solid line) and n = 1.0 x 10" 
cm~'^ (dashed line)^. (c) T2 as a function of tem- 
perature for GaSb with n = 2.0 x 10^^ cm^'^ (solid 
line) and n = 4.0 x 10^^ cm"'^ (dashed line), (d) 
H YS T for GaSb with n = 2.0 x 10" cm'^ (solid 
line) and n = 4.0 x 10^^ cm^"^ (dashed hne)^. 
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FIG. 3: Spin splitting of the lowest conduc- 
tion band for bulk III-V semiconductors at room 
temperature. (a) A£' vs k[iio] and (b) IS.E 
vs i^cxcoss for GaAs (solid line), GaSb (dashed 
line), InAs (dot-dashed line), and InSb (dot- 
ted Hne). (c) ^E vs kjno] and (d) IS.E vs 
-E-oxccss for GaP (solid line) and InP (dashed line). 
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FIG. 5: The fourteen-band electronic band 
structure near the P point of the Bril- 
liouin zone at 300K. (a) InAs. (b) InSb. 
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FIG. 4: Temperature dependence of mobility and 
coherence time in InAs and InSb. (a) T2 as 
a function of temperature for InAs with n ~ 
1.7 X IQi*' cm-3. (b) /x vs T for InAs^s. (c) 
T2 as a function of temperature for InSb with 

(d) ^ vs T for InSb^s. 



FIG. 6: Temperature dependence of mobility and 
coherence time in InP. (a) T2 as a function of tem- 
perature^" for n = 2.0 X 10^^ cm~^. (b) ^ vs T. 
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FIG. 7: The fourteen-band electronic 

band structure near the F point of 
the BriUiouin zone at 300K for InP. 



FIG. 9: The fourteen-band electronic band struc- 
ture of bulk II- VI semiconductors near the F point 
of the BriUiouin zone at BOOK, (a) CdSe. (b) ZnSe. 
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FIG. 8: Temperature dependence of mobility and 
coherence time in bulk II- VI semiconductors, (a) 
T2 as a function of temperature for CdSe with 
n = 4.0 X IQi^ cm-3. (b) /i vs T for CdSeSa. 
(c) T2 as a function of temperature for ZnSe with 

(d) ^ vs T for ZnSc™. 
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FIG. 10: Energy dependence of 1^^3(1, £^) of 
III-V semiconductors at room temperature for 
bulk GaAs (solid line), GaSb (long dashed 
line), InAs (long dot-dashed line), InSb (dot- 
ted line), InP (dot-dot dashed line), ZnSe (short 
dashed line), and CdSe (short dot-dashed line). 
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FIG. 11: Band-edge diagram (a) and band struc- 
ture (b) for a 75-A GaAs/100-A Alo.4Gao.6As 
(OOl)-oriented quantum well at 300 K. The la- 
bels c, I, h, and so identify the conduc- 
tion, light-hole, heavy-hole, and spin-orbit split- 
off hole band edges, respectively, at the zone 
center of the bulk constituent semiconductors. 
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FIG. 12: Spin splitting of the lowest conduction 
band for quantum wells and superlattices at room 
temperature, (a) AE vs krj^xo] ^n.d (b) AE vs -Eoxcess 
for a 75-A GaAs/100-A Alo.4Gao.6As QW (solid 
line), a thin-layer 21.2-A InAs/36.6-A GaSb SL 
(dashed Hue), an 70-A Ino.53Gao.47As/97-A InP QW 
(dot-dashed line), and a 80-A GaSbo.8i/80-A AlSb 
QW (dotted hne). (c) AE vs k[no] and (d) AE 
vs i^excess for a 150-A GaAs/100-A Alo.03Gao.97As 
QW (solid line), a 150-A GaAs/100-A Alo.4Gao.6As 
QW (dashed hne), a lOG-A Ino.13Gao.87As/2OO- 
A GaAs QW (dot-dashed), and a 95- A 

In0.57Ga0.43As0.93P0.07/75-Aln0.87Ga0.13As0.29P0. 71 
QW (dotted hne), and a 40-A Ino.20Gao.80As/lOO- 
A GaAs/40-A Ino.2oGao.8oAs/200- 

A GaAs QW (dot-dot-dashed hne). 
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FIG. 13: Ti as a function of (a) mobility fi, (b) 
average transport time Tp, and (c) carrier density 
n, for a 75-A GaAs/100-A Alo.4Gao.6As QW at 
room temperature. Closed circles represent the re- 
sults of experiments^^. The non-perturbative the- 
ory reults with OP (solid lines and open squares) 
and NI scattering (dashed lines and open trian- 
gles) are shown, as well as the DK theory re- 
sults (dot-dashed lines), (d) The measured car- 
rier density and mobility at room temperature^. 



FIG. 14: Electron spin coherence times and spin re- 
laxation rates as a function of electron confinement 
energy and quantum- well width at room temperature 
for (25-300)-A GaAs/100-A Alo.4Gao.6As quantum 
wellsi^ with electron density n = 2.686 x 10^^ cm~^ 
and electron mobility ^ = 800 cm^/Vs. Closed cir- 
cles represent the results of experiments. Solid and 
dashed lines respectively represent the results of the 
non-perturbative theory for OP and NI scattering. 
Dot-dashed lines represent results of the DK theory. 
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FIG. 15: Electron spin coherence times and spin 
relaxation rates as a function of electron confi- 
ment energy and quantum-well width at room tem- 
perature for a set of (25 - 300)-A GaAs/150-A 
A-lo.35Gao.65As quantum wells^^ with electron den- 



sity n 



1 X 10" cm 



Glosed circles rep- 



resent the results of experimental measurements, 
while solid and dashed lines represent the numer- 
ical calculations with /i = 4600 cm^/Vs for optical 
phonon and neutral impurity scattering, respectively. 



FIG. 16: Electron spin coherence times as a function 
of temperature. Glosed squares, circles, and triangles 
represent the experimental measurements for a 60- 
A GaAs/150-A Alo.35Gao.65As, a 100-A GaAs/150- 
A Alo.35Gao.65As, and a 200-A GaAs/150-A 
Alo.35Gao.65As quantum well, respectively, and the 
corrresponding calculated Ti's are respectively rep- 
resented by solid, dashed, and dot-dashed lines. 
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FIG. 17: Electron spin coherence times and spin re- 
laxation rates of (25 — 100)-A Ino.53Gao.47As/97-A 
InP QWs^^ as a function of electron confinement 
energy and quantum-well width at room temper- 
ature with ^ — 6700 cm^/Vs for non-degenerate 
(n = 6 X 10^^ cni^"^) and degenerate (n = 
5 X 10^* cm~^) electron densities. Closed cir- 
cles represent the results of experiments. Solid 
(dot-dot-dashed) and dashed (dot-dashed-dashcd) 
lines respectively represent the results of the non- 
perturbative theory for OP and NI scattering with 
non-degenerate (degenerate) electron density. Dot- 
dashed (dotted) lines represent results of the DK the- 
ory for non-degenerate (degenerate) electron density. 
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FIG. 18: Band-edge diagram (a) and band struc- 
ture (b) for a 70-A Ino.53Gao.47As/97-A InP 
(OOl)-oriented quantum well at 300 K. The la- 
bels c, I, h, and so identify the conduc- 
tion, light-hole, heavy-hole, and spin-orbit split- 
off hole band edges, respectively, at the zone 
center of the bulk constituent semiconductors. 
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FIG. 19: Band-edge diagram (a) and band struc- 
ture (b) for a 95-A Ino.57Gao.43Aso.93Po. 07/75- 
A Ino.87Gao.13Aso.29Po. 71 (OOl)-oriented quantum 
welii^ at 300 K. The labels c, I, h, and so iden- 
tify the conduction, light-hole, heavy-hole, and spin- 
orbit split-off hole band edges, respectively, at the 
zone center of the bulk constituent semiconductors. 
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FIG. 20: Ti as a function of confinement energy Ei^ 
at 300 K for (a) (30 - 100)-A GaSb/AlSb QWsii 
and (b) (30- 100)-A GaAsSb/AlSb QWe^^ with n = 
2 X 10^^ cm-3 and fi = 2000 cmVVs. Sohd and 
dashed lines respectively represent the results of the 
non-perturbative theory for OP and NI scattering. 
Dot-dashed lines represent results of the DK theory. 
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FIG. 21: Band-edge diagram (a) and band structure 
(b) for a 80-A GaSb /80-A AlSb (OOl)-oricnted quan- 
tum well at 295 K. The labels c, /, /i, and so iden- 
tify the conduction, light-hole, heavy-hole, and spin- 
orbit split-off hole band edges, respectively, at the 
zone center of the bulk constituent semiconductors. 



FIG. 22: Band-edge diagram (a) and band 
structure (b) for a 51-A Gao.igAso.siSb /80- 
A AlSb (OOl)-oriented quantum weU at 295 K. 
The labels c, I, h, and so identify the conduc- 
tion, light-hole, heavy-hole, and spin-orbit split- 
off hole band edges, respectively, at the zone 
center of the bulk constituent semiconductors. 
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FIG. 23: Electron spin coherence times. (a) 
Calculated Ti as a function of temperature T 
for a 80-A GaSb/300-A Alo.4Gao.6Sb QWB with 
NI scattering. (b) Measured mobility as a 

function of temperature-^^. Solid and dashed 
lines are for carrier density n = 2.11 x 10^^ 



cm 



and 



= 2.76 X 



respectively. 
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FIG. 24: Band-edge diagram (a) and band structure 
(b) for a 21.2-A InAs/36.6-A GaSb (OOl)-oriented 
strained-layer superlattice at 300 K. The labels c, I, 
h, and so identify the conduction, light-hole, heavy- 
hole, and spin-orbit split-off hole band edges, re- 
spectively, at the zone center of the bulk constituent 
semiconductors. The conduction, light-hole, heavy- 
hole, and spin-orbit split-off hole superlattice states 
are labelled by Ci, LHi, HHi, and SOi, respectively. 
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FIG. 25: Band-edge diagram (a) and band struc- 
ture (b) for a 75-A InAs/200-A AlSb (001)- 
oriented quantum well at 300 K. The labels c, 
/, and h identify the conduction, light-hole, and 
heavy- hole band edges, respectively, at the zone 
center of the bulk constituent semiconductors. 
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FIG. 26: T\ as a function of (a) confinement en- 
ergy _Bio and (b) quantum-well width at 300 K for 
(25 - 300)-A InAs/AlSb QWs with n = 2.00 x lO^^ 
cm^'^ and [i — 2000 cm-^/Vs. T\ as a function of (c) 
mobility with n = 2 x 10^^ cm^"^ and (d) electron 
density with [i = 2000 cm^/Vs at 300 K for a 75- 
A InAs/AlSb QW. Solid and dashed lines represent 
the results for OP and NI scattering, respectively. 
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27: Band-edge diagram (a) and band struc- 
(b) for a 100-A Ino.13Gao.87As/2OO-A GaAs 



(OOl)-oriented quantum wellSi at 300 K. The 
labels c, I, h, and so identify the conduc- 
tion, light-hole, heavy-hole, and spin-orbit split- 
off hole band edges, respectively, at the zone 
center of the bulk constituent semiconductors. 
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FIG. 28: Mobility and carrier-density dependence 
of spin coherence time in spin-LED structures at 
300K. (a) Ti vs ^JL with n = 1 x lO" cm^^ 
and (b) Ti vs n with /i = 2000 cm^/Vs for 
a 100-A Ino.13Gao.87As/2OO-A GaAs (OOl)-oriented 
quantum wellSl. (c) Ti vs /i with n = 1 x 
10" cm-3 and (d) Ti vs n with /i = 2000 
cmVVs for a 40-A Ino.20Gao.80As/lOO-A GaAs/40- 
A Ino.2oGao.8oAs/200-A GaAs (OOl)-oriented quan- 
tum well^. Solid and dashed lines respectively 
represent the results for OP and NI scattering. 
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FIG. 29: Band-edge diagram (a) and band 
structure (b) for a 40-A Ino.20Gao.80As/lOO- 
A GaAs/40-A Ino.2oGao.8oAs/200-A GaAs (001)- 
oriented quantum wcll^^ at 300 K. The labels 
c, Z, and h identify the conduction, light-hole, 
and heavy-hole band edges, respectively, at the 
zone center of the bulk constituent semiconductors. 



FIG. 30: Band-edge diagram (a) and band struc- 
ture (b) for a 150-A GaAs/100-A Alo.03Gao.97As 



(OOl)-oriented quantum welisa at 300 K. The 
labels c, ?, /i, and so identify the conduc- 
tion, light-hole, heavy-hole, and spin-orbit split- 
off hole band edges, respectively, at the zone 
center of the bulk constituent semiconductors. 
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FIG. 31: Band-edge diagram (a) and band struc- 
ture (b) for a 150-A GaAs/100-A Alo.4Gao.6As 
(OOl)-oricnted quantum well at 300 K. The la- 
bels c, I, h, and so identify the conduc- 
tion, light-hole, heavy-hole, and spin-orbit split- 
off hole band edges, respectively, at the zone 
center of the bulk constituent semiconductors. 
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FIG. 32: Mobility and carrier-density dependence of 
spin coherence time in spin-LED structures at room 
temperature. Ti vs ^ with n = 1 x 10^^ cni~^ 
(a) and Ti vs n with ^ = 2000 cmVVs (b) for 
a 150-A GaAs/lOQ-A Alo.03Gao.97As (OOl)-oriented 
quantum well^. Ti vs /i with n — 1 y. 10^^ cm~^ 
(c) and Ti vs n with ^ = 2000 cmVVs (d) for 
a 150-A GaAs/100-A Alo.4Gao.6As (OOl)-oriented 
quantum well. Solid and dashed lines respectively 
represent the results for OP and NI scattering. 
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FIG. 33: Calculated T\ as a function of (a) mobil- 
ity //, (b) average transport time Tp, and (c) Tem- 
perature T, for a 50-A GaAs/100-A Gao.51Ino.49P 
QW2i with carrier density n = 2.67 x lO" 
cm~'^. (d) Measured mobility as a function of 
temperature^,. Solid and dashed lines respectively 
represent the results for OP and NI scattering. 
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FIG. 34: Calculated T\ as a function of (a) mo- 
bility ^, (b) average transport time Tp, and (c) 
Temperature T, for a 105-A CdSe/500-A ZnSe 
QW^ with carrier density n — 7.93 x 10^^ cm"'^. 
(d) Extrapolated mobility as a function of tem- 
perature2^. Solid and dashed lines respectively 
represent the results for OP and NI scattering. 



45 





;(a) 




- 


3.20 


- 




- 


3.00 






- 














r- 


- 






C 




2.80 


_ 


N 


_ 




- u 




CD - 




- C/5 


cz 


&n - 




- a 


■a 


C - 




-N 


U 


N - 


2.60 














- 


> 




C - 


0) 




(; 




>C2.40 


^ 




J 


^ 




h 


_ 


t^'O.OO 


- 




- 






I 


h : 






/ - 


-0.25 






- 






so 


- 


-0.50 


- 




so - 




1 


1 , 


1 1 1 1 1 1 



-150 150 

Position (A) 




k||[110](r) 



FIG. 35: Band-edge diagram (a) and band struc- 
ture (b) for a 105-A Cdo.24Zno.76Se/5OO-A ZnSe 
(OOl)-oriented quantum weliii at 300 K. The 
labels c, I, h, and so identify the conduc- 
tion, light-hole, heavy-hole, and spin-orbit split- 
off hole band edges, respectively, at the zone 
center of the bulk constituent semiconductors. 
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FIG. 36: Calculated Ti as a function of (a) mo- 
bility ^, (b) average transport time Tp, and (c) 
Temperature T, for a 105-A Cdo.24Zno.76Se/5OO-A 
ZnSe QW^ with carrier density n = 7.44 x IQi^ 
cm~^. (d) Extrapolated mobility as a function 
of temperature^^. Solid and dashed lines respec- 
tively represent the results for OP and NI scattering. 
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FIG. 37: Energy dependence of electron spin 
precession vector for III-V quantum wells and su- 
perlattices at room temperature, (a) i^'^^{l,E) and 
(b) ^23(1,^) for a 75-A GaAs/100-A Alo.4Gao.6As 
QW (solid line), a thin-layer 21.2-A InAs/36.6- 
A GaSb SL (dashed line), a thin-layer 75-A 
InAs/200-A AlSb QW (dot-dot dashed hne), an 
70-A Ino.53Gao.47As/97-A InP QW (dot-dashed 
line), and a 80- A GaSbo.8i/80-A AlSb QW (dotted 
line). The short-dashed line is the DK approx- 
imation for the 75-A GaAs/100-A Alo.4Gao.6As 
quatum well. (c) hli{l,E) and (d) ^2^3(1,^) 
for a 150-A GaAs/100-A Alo.03Gao.97As QW 
(solid line), a 150-A GaAs/100-A Alo.4Gao.6As 
QW (dashed line), a 100-A Ino.13Gao.87As/2OO- 
A GaAs QW (dot-dashed), and a 95- A 

In0.57Ga0.43As0.93P0.07/75-Aln0.87Ga0.13As0.29P0. 71 
QW (dotted hne), and a 40-A Ino.20Gao.80As/lOO- 
A GaAs/40-A Ino.2oGao.8oAs/200- 

A GaAs QW (dot-dot-dashed line). 



